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Chapter  I 
INTRODUCTION 

A cushioning  system  is  an  interface  device  which  isolates  a pro- 
tected item  from  the  shock  loads  which  develop  in  its  external  environ- 
ment. A cushioning  system  is  required  when  the  protected  item  cannot 
survive  the  shock  loads  imposed  by  the  external  environment,  unless  these 
loads  are  attenuated  by  some  form  of  cushioning  system.  The  design  of 
a cushioning  system,  therefore,  involves  stipulations  concerning  both 
the  item  to  be  protected  and  the  external  environment.  These  two  stipu- 
lations of  how  much  shock  an  item  can  survive  and  how  much  shock  it  is 
expected  to  experience  in  its  environment  are  questions  that  must  be 
answered  in  the  design  of  any  cushioning  system,  commercial  or  military. 
It  is  the  concern  for  extreme  external  environments  that  is  emphasized 
most  in  the  military  designs. 

The  Military  Environment 

A cushioning  system  designed  for  military  use,  such  as  a shock 
isolating  component  of  a missile  or  missile  system,  is  required  to  per- 
form its  basic  functions  throughout  its  operating  life.  Current  mili- 
tary policy,  uased  on  the  need  for  worldwide  deployment,  dictates  that 
military  materiel  be  capable  of  withstanding  the  rigors  of  deployment 
on  a worldwide  basis.  Consequently,  military  materiel,  including 
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cushioning  systems,  must  withstand  all  the  environments  to  which  they 
are  exposed  and  continue  to  operate  satisfactorily  if  the  materiel  is 
to  be  considered  acceptable. 

The  environment  in  which  military  materiel  operates  is  frequently 
quite  severe  and  may  produce  substantial  damage.  Military  operations 
encompass  all  the  geographic  regions  and  the  peculiar  aspects  of  their 
environments.  Each  environment  introduces  hazards  peculiar  to  geo- 
graphic zones  which  material  will  encounter  when  deployed. 

The  tropics  present  a hot,  humid  atmosphere,  subjectin''  materiel 
to  what  is  commonly  referred  to  as  "jungle-rot".  Excessive  rainfall, 
high  humidity,  heat,  and  fungus,  in  addition  to  the  prevalence  of  vermin, 
insects,  and  reptiles  present  problems  relating  to  materiel  protection 
and  to  the  safety  of  personnel. 

The  arctic  regions  subject  materiel  to  a very  cold  atmosphere. 
Temperatures  of  -65°F  are  common  and  temperatures  of  -85°F  have  been 
recorded  in  underground  ammunition  storage  "igloos".  The  physical 
characteristics  of  components,  in  particular  resilient  cushioning  sys- 
tems, are  radically  altered  under  these  conditions  and  their  operational 
integrity  is  jeopardized. 

Also,  high-speed,  fast-climbing  aircraft  introduce  rapid  and 
severe  variations  in  temperature  and  pressure,  and  provisions  must  be 
made  to  withstand  the  effects  of  these  changes,  when  required. 

Cushioning  Systems  in  the  Military 

In  evaluating  the  environmental  factors,  the  military  has  estab- 
lished certain  extreme  limits  on  the  particular  conditions  that  have 
become  accepted  as  a basis  for  worldwide  environments.  Temperature 
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extremes,  which  produce  dramatic  changes  in  the  resiliency  of  cushion- 
ing materials,  are  important  in  cushion  system  design.  The  world- 
wide temperature  extremes  are  generally  defined  as  -65°F  and  160°F; 
cushioning  systems  in  military  applications  are  expected  to  perform 
their  shock  interface  function  under  these  adverse  environments. 

Many  of  the  military  cushioning  systems  use  some  form  of  bulk 
cushioning  material  (e.g.  plastic  foam)  as  the  cushioning  system.  This 
provides  a low-cost,  lightweight  cushioning  system  that  is  easily  incor- 
porated into  the  design;  however,  designers  have  some  reservations  con- 
cerning the  use  of  bulk  cushioning  systems  in  the  military.  One  limit- 
ing factor  is  the  difficulty  the  designer  has  in  predicting  the  response 
of  bulk  cushioning  materials  when  subjected  to  the  extreme  environmental 
conditions  encountered  under  military  deployment.  The  temperature 
sensitivity  of  these  cushioning  materials,  which  causes  variations  in 
the  impact  response,  is  also  a primary  concern.  It  was  shown  in  a 
recent  study  of  thermoplastic  foam  cushioning  systems  that  temperature 
had  a significant  effect  on  cushion  system  response  [1],  Therefore, 
temperature  effects  must  be  factored  into  the  design  of  cushion  systems 
using  bulk  cushioning  materials. 

To  fully  utilize  bulk  cushioning  systems,  designers  require  suffi- 
cient information  to  accurately  predict  response  variations.  The  cur- 
rent practice  is  to  provide  the  designer  with  cushioning  data  for  each 
type  and  thickness  of  cushioning  material.  These  data  are  provided  in 
the  form  of  dynamic  cushioning  curves  as  prescribed  in  the  current 
theory  of  cushioning  design  (Chapter  II) . For  any  particular  shock 
isolation  system  design  program,  the  designer  is  generally  given  a 
maximum  allowable  fragility  level  which  the  protected  item  is  permitted 


to  experience.  Also  the  particular  organization  involved  will  have  an 
established  testing  policy  defining  appropriate  impact  tests.  These 
parameters  provide  the  basis  for  the  design  of  the  shock  mitigation 
system.  Given  a satisfactory  prediction  of  impact  response  for  various 
candidate  cushioning  materials,  a designer  can  use  this  prediction  to 
select  the  appropriate  cushioning  scheme  to  meet  the  particular  design 
requirements. 


Research  Objective 

The  research  objective  of  this  study  is  to  develop  a reliable 
impact  response  model  that  accurately  predicts  the  dynamic  cushion- 
ing performance  of  bulk  cushioning  systems.  A secondary  objective  is 
to  develop  an  optimization  technique  that  utilizes  the  model  in  deter- 
mining an  optimal  cushioning  system  design. 

In  answering  the  research  objective  a systematic  study  of  back- 
ground material  was  conducted.  The  current  theory  of  cushioning  design 
is  discussed  in  Chapter  TI.  The  ingredients  of  the  General  Model  of 
impact  response  are  identified  in  Chapter  III  and  the  basis  of  the 
underlying  structure  of  the  developed  model  is  based  upon  viscoelastic 
theory. 

The  modeling  process  and  the  analysis  techniques  used  to  formulate 

the  General  Model  of  impact  response  are  also  discussed  in  Chapter  III. 

Chapter IV  presents  the  two  finalized  models,  a General  Model  of  impact 

response  and  the  Minicel  Model,  the  model  of  impact  response  of  a 
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particular  cushioning  material  (Minicel  is  a 2 lb/ ft  cross-linked 
polyethylene  foam  manufactured  by  Hercules,  Inc). 
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The  validity  of  the  models  are  demonstrated  in  Chapter  V through 
a systematic  series  of  tests  and  analysis.  Finally,  in  answer  to  the 
secondary  objective  of  the  research,  an  optimization  technique  is  pre- 
sented in  Chapter  VI  that  generates  the  optimal  bulk  cushion  design. 

Th(»  optimization  technique  utilizes  the  predictive  capability  of  the 
General  Model  of  impact  response  to  determine  optimality  and  provide, 
as  output,  a set  of  dynamic  cushioning  curves  at  the  optimal  conditions. 


Chapter  II 


CUSHIONING  DESIGN  THEORY 

Anything  that  is  subject  to  movement  is  subject  to  mechanical 
damage  due  to  shock.  Shock:  may  be  defined  as  a sudden  change  in  direc- 
tion or  velocity  of  the  motion  of  a body.  The  magnitude  or  intensity 
of  shock  is  expressed  in  G's,  which  is  defined  in  terms  of  the  time 
rate  of  change  of  the  velocity  (acceleration),  and  is  measured  in 
feet/ second/ second.  Mathematically*  G's  = a/g  where  a is  the  accelera- 
tion experienced  by  the  body  and  g i>  the  acceleration  due  to  gravity 
(32.2  ft/sec/ sec). 

A given  body  in  a static  condition  has  one  gravity  unit  of  G act- 
ing on  it  and,  therefore,  exerts  one  G upon  its  support.  If  this  body 
is  raised  and  allowed  to  fcH  freely,  it  will  accelerate  in  its  fall, 
due  to  the  force  of  gravity,  until  it  collides  with  its  support  or  the 
earth.  The  stop  causes  the  body  to  experience  a sudden  deceleration 
that  can  be  expressed  in  terms  of  G's.  If  the  body  experiences,  upon 
impact,  a deceleration  of  20  times  that  of  its  static  condition,  it  is 
said  to  have  experienced  20  G's  of  deceleration.  A jet  pilot  experi- 
ences such  a condition  when  he  pulls  out  of  a dive.  He  must  not  exceed 
his  G limit  or  he  will  black  out,  and  if  the  aircraft  is  not  designed 
and  stressed  to  withstand  high  G's  during  such  maneuvers,  severe  damage 
will  occur  to  the  aircraft  and  it  may  crash.  This  same  situation  exists 
in  regard  to  fragile  objects.  The  fragility  level  of  an  object,  measured 
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in  G's,  is  its  ability  to  withstand  deceleration.  The  purpose  of  a 
cushioning  system  is  to  reduce  the  shock  transferred  to  the  protected 
object  to  a degree  below  its  fragility  level,  thus  protecting  it  against 
physical  damage. 


Cushioning  System  Design 

To  understand  how  a cushioning  material  functions  during  shock 
transfer,  one  can  consider  what  happens  when  a body,  is  dropped  onto  a 
rigid  surface  such  as  a concrete  floor.  At  impact  the  body  is  falling 
at  a velocity  (V  = Jl  gh  ) , where  V is  the  velocity  at  impact  in 
feet/second,  g is  the  acceleration  due  to  gravity  of  32.2  feet/second/ 
second,  and  h is  the  height  of  drop  measured  in  feet.  In  a very  short 
time  after  impact,  this  velocity  is  reduced  to  zero  in  a very  small 
distance.  Thus,  a rapid  decrease  in  velocity  occurs,  due  to  impact, 
and  the  body  is  subjected  to  a very  high  deceleration. 

If  the  same  situation  is  considered  except  that  the  body  is  drop- 
ped onto  a resilient  cushioning  material  which  rests  on  the  same  rigid 
surface,  the  body  has  the  same  velocity  at  impact.  However,  the  time 
required  for  the  velocity  to  be  reduced  to  zero  is  much  greater  than 
in  the  previous  situation;  the  rate  at  which  the  velocity  decreases  is 
considerably  less;  and  the  distance  traveled  afte  initial  contact  with 
the  cushion  until  the  time  the  velocity  is  reduced  to  zero  is  consider- 
ably greater.  Compared  to  the  first  situation,  the  body  is  subjected 
to  lower  G's.  The  resilient  cushioning  material  has,  therefore,  atten- 
uated the  shock  pulse  by  dissipating  the  kinetic  energy  present  in  the 
body  at  the  time  of  impact  over  a longer  time  period.  This  has  been 
expressed  mathematically  [2]  as  follows: 

6-  f\F“" 
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where 


G * acceleration  G - level 
t = shock  pulse  rise  tine  in  milliseconds 
h * drop  height  in  inches. 

Equation  (11-1)  shows  that  as  the  shock  pulse  rise  time  is  increased, 
the  G’s  experienced  by  the  body  are  proportionally  decreased. 

Increases  in  shock  pulse  rise  time  can  usually  be  accomplished  by 
increasing  the  thickness  of  the  resilient  cushioning  material.  For 
example,  the  G's  experienced  by  a body  falling  on  a tangentially  elastic 
cushioning  material  can  be  predicted  as  follows  [ 2] 


G = 


3.9  h 
T 


(II-2) 


where  T = thickness  of  cushion  in  inches. 

Equation  (II-2)  shows  that  the  predicted  G levels  are  reduced 
proportionally  with  increases  in  the  thickness  of  the  cushioning 
material.  If  the  cushion  thickness  is  increased,  then,  during  impact, 
the  excursion  envelope  of  the  body  is  increased  proportionally  and  the 
body  moves  through  an  additional  amount  of  cushioning  material  before 
it  comes  to  rest.  This  increase  in  excursion  directly  increases  the 
shock  pulse  rise  time  with  an  accompanying  reduction  in  G levels.  This 
reduction  is  consistent  wich  the  change  in  G levels  that  would  be  pre- 
dicted by  Equation  (II-l)  with  an  increase  in  shock  pulse  rise  time. 


History  of  Cushion  System  Design 

Equation  (II-2)  is  based  on  Mindlin's  work  in  1945  [3]  that  marked 
the  beginning  of  the  scientific  approach  to  cushion  design.  Mendlin's 
paper  was  followed  by  a number  of  discussions  [4-8]  by  author's  who 
adopted  his  procedures.  The  next  substantial  step  forward  was  the 


development  of  optimum  efficiency  design  points  by  Janssen  [ 9]  . Janssen 
utilized  material  properties  determined  by  quasistatic  means,  in  partic- 
ular, stress-strain  curves  determined  on  a conventional  compression 
tester  with  the  speed  of  compression  quite  slow  (not  more  than  2 inches/ 
minute).  He  derived  a cushion  factor,  "J" , for  several  cushioning 
materials  that  was  the  ratio  of  optimum  stress  to  optimum  strain.  Then 
G-level  can  be  predicted  on  the  basis  of 

G - J £ (II-3) 

where  J = Janssen's  cushion  factor. 

This  was  a significant  improvement  over  Equation  (II-2)  in  that 
the  J values  allowed  for  different  performance  factors  for  different 
materials.  However,  it  soon  became  apparent  that  single  curves  based 
on  the  static  stress-strain  characteristics  of  a material  did  not 
describe  the  dynamic  behavior.  Several  methods  of  presentation,  all 
involving  families  of  curves,  were  attempted.  Gradually,  Kerstner's 
approach  [10]  became  the  most  popular.  In  this  approach,  a family  of 
curves  is  drawn  for  each  material  thickness  at  each  drop  height,  plot- 
ting peak  dynamic  stress  (or  acceleration)  against  the  original  static 
stress.  A typical  set  of  such  curves,  taken  from  Humbert  and  Hanlon 
[ 11]  , is  shown  in  Figure  1. 

These  curves,  referred  to  as  dynamic  cushioning  curves,  are  gener- 
ated for  a particular  type  and  thickness  of  cushion  by  performing  drop 
tests  using  standard  weight  specimens  that  are  dropped  onto  the  cushion. 

The  static  stress  (og)  is  determined  by: 


(II-4) 


where  a is  the  static  stress  (psi),  W is  the  specimen  weight,  and  A is 
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the  footprint  of  the  specimen  in  the  cushion  (in.  ).  A different  curve  is 
required  for  each  drop  height,  thickness,  and  type  of  material.  The 
curves  are  a good  indication  of  the  protection  to  be  expected  for  a 
particular  cushion  scheme. 

The  current  practice  utilizes  the  Kerstner  procedure  in  that 
manufacturers  of  cushioning  products  demonstrate  the  cushioning  perfor- 
mance capabilities  of  their  particular  products  by  generating  a series 
of  dynamic  cushioning  curves  at  various  drop  heights  and  thicknesses  of 
cushions  (Figure  1).  Also,  Military  Handbook  -304,  "Cushioning  Design 
Handbook",  was  published  in  1964  which  provides  families  of  dynamic 
cushioning  curves  for  a large  number  of  frequently  utilized  cushioning 
materials.  Additional  work  was  done  by  Mustin  [12]  wno  used  a single 
value  for  correlating  each  dynamic  cushioning  curve  in  a family  of  curves 
such  as  Figure  1.  This  proved  to  be  an  over  simplification  in  develop- 
ing a general  model  of  impact  response,  as  will  be  seen  in  Chapter  IV. 

Reservations  Concerning  Current  Cushioning  System  Design 

In  recent  years,  equipment  designers  have  become  increasingly  aware 
of  the  detrimental  effects  of  extreme  temperature  upon  equipment  per- 
formance. Consequently,  there  are  now  included  in  the  qualification 
tests  of  equipment,  some  tests  conducted  at  temperatures  that  are  repre- 
sentative of  the  temperature  extremes  that  the  equipment  is  likely  to 
encounter . 

This  extreme  temperature  testing  has  received  substantial  atten- 
tion in  the  military,  where  the  range  of  temperatures  encountered  is 
quite  extreme,  and  the  failures  can  produce  catastrophic  results.  The 
transportation  of  military  equipment  is  one  area  of  concern  and  a study 
was  made  of  several  military  containers  that  used  bulk  cushioning  sys- 
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terns.  It  was  found  that  temperature  appeared  to  have  a significant 
effect  on  Impact  response  [ 1] . The  cushioning  systems  in  the  containers 
did  not  perform  properly  due  to  changes  in  the  performance  of  the  bulk 
cushioning  materials  that  were  induced  by  extreme  temperatures.  Conse- 
quently, the  items  packaged  in  the  containers  (guided  missile  systems 
and  system  components)  did  not  receive  adequate  protection,  and  the 
missile  system  reliability  was  compromised.  This  type  of  failure  is  a 
potential  problem  that  can  occur  in  very  expensive  equipment  and  pro- 
duce a malfunction  in  weapon  systems  that  compromises  the  combat  power 
of  a military  organization.  Also,  if  proper  failsafe  provisions  are 
not  incorporated  into  the  protection  of  ordinance  items,  the  safety  of 
any  of  the  personnel  that  handle  the  equipment  within  the  logistics 
system  is  jeopardized. 

Rather  than  incorporate  additional  protection  into  equipment,  it 
is  much  more  cost  effective  to  improve  the  reliability  of  cushioning 
systems.  This  can  be  done  if  a reliable  method  of  predicting  cushion- 
ing performance  can  be  developed. 

A technique  has  been  suggested  for  modifying  the  conventional 
dynamic  cushioning  curves  in  such  a manner  as  to  address  the  effect  of 
temperature  on  the  impact  response  of  bulk  cushioning  systems  [1].  The 
technique  utilizes  superimposed  dynamic  cushioning  curves  that  are  a 
super-positioning  of  the  dynamic  cushioning  curves  at  temperature 
extremes  upon  the  ambient  dynamic  cushioning  curve.  One  such  curve  is 
presented  in  Figure  2.  This  type  of  curve  demonstrates  the  effect  of 
temperature  for  a selected  set  of  conditions  and  provides  the  cushioning 
system  designer  with  the  capability  of  designing  cushioning  systems  with  a 
reduced  likelihood  of  design  failure  under  extreme  temperature  conditions. 
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Figure 


Chapter  ill 

BASIS  FOR  THE  IMPACT  REbrurtSE  liODEL 


A valid  model  of  impact  response  must  incorporate  all  parameters 
that  are  expected  to  have  a significant  effect  on  impact  response. 

Temperature  has  a significant  effect  on  impact  response  [1],  and  it  is 
postulated  that  viscoelastic  theory  can  be  utilized  to  formulate  a 
model  of  impact  response  that  incorporates  temperature  effects.  The 
current  design  practice  for  predicting  impact  response  is  predicated  on 
dynamic  cushioning  curves.  Dynamic  cushioni*  3 curves  do  not  account 
for  temperature  effects  on  impact  response.  To  improve  the  predicta- 
bility of  cushioning  systems,  a model  of  impact  response  must  account 
for  temperature  effects.  The  temperature  effects  would  be  expected  to 
be  the  most  dramatic  as  the  temperature  tends  towards  the  extremes; 
therefore,  the  temperature  extremes  encountered  by  a cushioning  system 
are  of  particular  concern  in  modeling  impact  response. 

Temperature  Extremes  in  the  Model 

Army  Regulation  70-38,  "Research,  Development,  Test  and  Evaluation 
of  Materiel  for  Extreme  Climatic  Conditions,"  requires  that  the  extreme 
external  environments  that  are  likely  to  be  encountered,  be  considered 
in  the  design  of  Army  materiel.  The  temperature  extremes  to  be  used 
in  the  design  are  defined  in  the  AR  according  to  the  intended  deploy- 
ment of  the  material  being  designed.  Figures  3 and  4 (from  AR  70-38) 
present  maps  of  the  extreme  temperature  conditions  to  be  used  in  the 


Figure  4.  Distribution  of  absolute  maximum  temperatures. 
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design  of  Army  materiel.  For  example,  for  worldwide  deployment,  the 
operational  temperature  range  is  given  as  -65°  to  160°F.  However,  for 
continental  United  States  deployment  only,  the  range  is  from  -30°  to 
145°F.  Army  materiel  is  expected  to  function  properly  within  these 
ranges  and  consequently  the  qualification  testing  of  Army  materiel  is 
conducted  at  these  extremes. 

A research  effort  was  initiated  at  the  Army  Missile  Command  (MICOM) 
in  1973  to  develop  superimposed  dynamic  cushioning  cqrves  that  address 
the  effect  of  temperature  on  the  impact  response  of  cushioning  materials. 

The  drop  tests  were  conducted  at  MICOM  and  the  results  analyzed  by 
the  University  of  Alabama  in  Huntsville  (UAH)  under  a supporting  research 
contract  [ 13] . The  initial  experimentation  was  conducted  on  Minicel 

3 

material,  a cross-linked  polyethylene  foam  material  with  a 2 pound/ foot 
density  manufactured  by  Hercules,  Inc. 

The  drop  test  program  that  was  conducted  on  the  Minicel  material 
used  drop  heights  of  12,  18,  24,  and  30  inches;  temperatures  of  -65°, 

70°,  and  160°F;  cushion  thicknesses  of  1,  2,  and  3 inches;  and  static 
stress  levels  that  varied  from  0.04  to  5.0  psi.  The  G-level  response 
and  shock  pulse  duration  were  recorded  for  each  of  2736  drop  tests.  An 
automated  data  handling  system  that  included  outlier  tests  and  other 
statistical  analyses  was  developed  and  used  to  analyze  the  Minicel  data. 

This  experimental  effort  resulted  in  the  generation  of  a data  base 
of  2409  statistically  valid  data  points.  A family  of  second  order 
polynomial  equations  was  found  to  be  the  best  predictor  of  impact 
response  of  the  Minicel  cushioning  material.  The  data  are  given  in 
Appendix  A together  with  tl  e families  of  regression  polynomials  for  the 
12  inch  and  24  inch  drop  heights,  and  the  correlation  coefficients. 


Theoretical  Basis  of  the  Model 


The  construction  of  an  impact  response  model  for  cushioning  mate- 
rials at  varying  temperatures  requires  the  development  of  a functional 
relationship  of  the  variables.  The  required  relationshio  can  be  expres- 
sed mathematically  as  follows: 

G = F(os,  l,  e,  h)  (III-l) 

where 

G = acceleration  G-level 

a = static  stress  in  psi 
s r 

T = thickness  of  cushion  inches 
6 = cushion  temperature  in  °F 
h » drop  height  in  inches. 

After  considerable  research,  it  was  determined  that  the  theory  of 
viscoelasticity  could  be  utilized  to  provide  a theoretical  basis  for 
the  model.  Viscoelastic  theory  recognizes  cushioning  materials  as 
belonging  to  a class  of  materials  which  have  mechanical  properties  that 
are  common  to  perfect  solids  and  perfect  liquids.  Various  theories 
have  been  developed  over  the  past  century  for  describing  the  behavior 
of  perfect  solids  and  perfect  liquids.  Among  these,  the  oldest  theories 
are  the  classical  theory  of  elasticity  and  the  theory  of  hydrodynamics. 
The  classical  theory  of  elasticity  deals  with  the  behavior  of  solids 
for  which  the  stress  is  directly  proportional  to  the  strain  but  is 
independent  of  rate  of  strain.  This  type  of  solid  is  known  as  a Hookian 
solid  (perfect  elastic  solid).  The  theory  of  hydrodynamics  describes 
the  behavior  of  perfect  viscous  liquids  for  which,  in  accordance  with 
Newton's  viscosity  law,  the  stress  is  directly  proportional  to  the  rate 
of  strain,  rather  than  the  strain  itself.  In  certain  instances,  solids 


and  liquids  may  have  their  stress  related  to  strain,  rate  of  strain,  and 
higher  time  derivatives  of  strain.  Behavior  of  such  materials  is  termed 
viscoelastic  when  the  stress  is  linearly  proportional  to  strain.  Mate- 
rials whose  behavior  is  viscoelastic  display  solid-like  and  liquid-like 
characteristics  [ 14] . 

Mathematical  models  have  been  formulated  for  viscoelastic  materials 
that  have  validity  over  both  short  and  long  time  periods;  for  example, 
Mustin  [ 12]  gives  the  creep  and  relaxation  functions  (long  term  behavior) 
for  a number  of  simple  mathematical  models  made  up  of  simple  spring 
elements  that  are  assumed  linear  and  massless,  and  of  dashpots  in  which 
a piston  is  moving  through  a liquid  that  obeys  Newton's  law  of  viscosity 
(velocity  is  proportional  to  strain).  The  short  term  stress  law  and  the 
long  term  creep  and  relaxation  functions  are  given  in  Table  I.  Creep 
TABLE  I.  SOME  CREEP  AND  RELAXATION  FUNCTIONS 


Function 

Linear 

Dashpot 

Voigt  Solid 

Maxwell  Solid 

Pictorial 

E 

c 

—> 

d€ 

a - c- 

_e  nc 

representation 
Stress  law 

a = Ee 

0 " Ee  + c dt 

£ + .1  do  m de 
c E dt  dt 

Creep  function 

i • 

“ • H(t) 

I (1  - e'Et/c)  . H(t) 

(!+  ;)■"« 

Relaxation 

function,  R(t) 

E • H(t) 

c • 6(t) 

E • H(t)  + c 5(t) 

(Ee‘Et/c)  • H(t) 

Notes: 


E = spring  constant 
o = stress 
e = strain 

c = damping  coefficient 
H(t)  = unit  step  function 
S(t)  = Dirac  delta  function 
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functions  are  shown  in  Figure  5 while  Figure  6 illustrates  the  relaxa- 
tion behavior.  The  appearance  of  the  relaxation  functions  indicates 
that  the  stress  is  infinite  at  the  instant  that  strain  occurs.  This  is 
due  to  the  dashpot  element  which,  unlike  a spring,  cannot  give  a finite 
instantaneous  strain  response  to  a finite  instantaneous  force  change. 
Therefore,  an  infinite  force  is  required  to  produce  a finite  instan- 
taneous strain. 


(a) 


Z' 

l 

_J 

TIME 

(b) 


Figure  5.  Behavior  of  some  simple  creep  functions: 
(a)  Maxwell  solid,  and  (b)  Voigt  solid. 


Figure  6.  Behavior  of  some  simple  relaxation  functions  (R(t)): 
(a)  Maxwell  solid,  and  (b)  Voigt  solid. 


These  relationships  are  expressed  compositely  in  constitutive 
equations  which,  in  viscoelastic  theory,  describe  the  response  of 
materials  to  mechanical  excitation.  The  constitutive  equation  in  con- 
junction with  the  energy,  momentum,  and  continuity  equations  permits 
the  prediction  of  complicated  responses  to  complicated  excitations, 
like  the  flow  of  liquids  in  irregularly  shaped  channels  or  the  deflec- 
tions of  beams  of  varying  cross  sections  under  complicated  loading 


conditions . 
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There  are  two  approaches  to  obtaining  the  constitutive  equation 
of  a material:  it  can  be  obtained  either  experimentally  with  the  help 

of  some  simple,  well-defined  tests  (like  the  stress-strain  curve)  or 
the  response  of  the  material  to  some  strain  or  stress  history  can  be 
calculated  with  the  help  of  a model  describing  its  structure.  Hooke's 
law  is  an  example  for  the  first  group,  the  phenomenological  equations. 
It  is  based  solely  on  experimental  observation.  The  best  known  example 
in  the  latter  category  of  structural  theories  is  the  kinetic  theory  of 
rubber  elasticity  predicting  the  elastic  response  of  vulcanized  elasto- 
mers from  their  structure. 

While  the  phenomenological  constitutive  equations  describe  the 
results  of  experimentally  obtained  data  and  are  usually  applied  to  pre- 
dict more  complicated  behavior,  the  structural  ones  offer  an  insight 
into  structure-property  relations  [ 15] . 

Since  temperature  effects  are  required  in  the  constitutive  equa- 
tion that  models  impact  response,  one  must  consider  that  portion  of 
viscoelastic  theory  which  accounts  for  the  behavior  of  materials  at 
varying  temperatures.  Most  of  this  theory  arises  from  the  considera- 
tion of  the  behavior  of  the  molecule.  The  molecule  may  be  visualized 
as  a long  curled  elastic  chain  which  may  have  cross-linking  with  other 
molecular  chains.  Since  molecular  activity  is  a function  of  tempera- 
ture, strain  response  (the  summation  of  individual  molecule  responses) 
is  also  a function  of  temperature.  On  this  basis,  in  the  glassy  zone 
(the  zone  below  the  temperature  at  which  the  polymer  structure  starts 
to  become  brittle),  molecular  "curling  and  uncurling"  cannot  occur 
rapidly  enough  to  follow  the  stress.  In  the  rubbery  zone  (the  zone 
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above  the  temperature  at  which  the  polymer  structure  becomes  rubbery) , 
curling  and  uncurling  are  in  phase  with  the  stress  which  is  not  condu- 
cive to  energy  dissipation. 

Let  the  long  term  elasticity  of  a material  at  a given  temperature,  , be 
plotted  against  the  natural  logarithm  of  time.  Suppose,  now,  that  the 
curve  shifts  horizontally  to  the  right  as  the  temperature  is  lowered. 
This  situation  is  illustrated  in  Figure  7 for  three  temperatures. 


Figure  7.  Shift  in  elasticity  characteristics 
as  a function  of  time  and  temperature. 

The  behavior  sketched  in  Figure  7 is  typical  of  many  materials  utilized  as 
cushions.  Materials  which  behave  in  this  manner  are  called  thermorheo- 
logically  simple.  Amorphous  polymers  behave  in  this  manner  while  poly- 
mers with  crystalline  structures  do  not.  Most  cushioning  materials  may 
be  considered  amorphous.  Teflon  and  plasticized  polyvinyl  chloride  are 
not  amorphous,  but  they  are  rarely  used  as  cushions  [16], 

Viscoelastic  Considerations 

A constitutive  equation  that  relates  the  maximum  G-level  experi- 
enced by  a body  when  subjected  to  a deceleration  impact  into  a cushion- 
ing material  is  required  for  the  impact  response  model  stated  in  Equa- 
tion (III-l) . A partial  basis  is  found  in  a paper  on  the  viscoelastic 


properties  of  thermorheologically  simple  cushioning  materials  by  Cost 
[ 17] . Cost  considers  a viscoelastic  Kelvin  body  which  is  identical  to 
the  Voight  model  shown  in  Table  1.  The  stress  laws  are  identical  with 
Figure  7 except  for  notation,  in  that  Cost  uses  K and  t)  for  E and  c, 
respectively. 

Consider  a body  of  mass  M falling  under  the  influence  of  gravity 

from  a height  h and  impacting  a cushion  material  of  area  Ac  and 

thickness  T.  The  static  stress  a is  defined  as  the  weight  of  the  body 

s 

divided  by  the  area  of  the  cushion  A . 

c 

Assume  the  cushion  material  behaves  as  a viscoelastic  Kelvin  body 
whose  stress  strain  relation  can  be  expressed  as 

a = Ke  + Tie  (III-2) 

where  q and  K are  material  properties.  When  Equation  (III-2)  is  applied 

to  the  problem  under  consideration,  0 becomes  the  static  stress  in  a 

bulk  cushion,  o , and  e,  the  unit  strain,  can  be  expressed  in  terms  of 
s 

the  displacement  of  the  upper  surface  of  the  cushion  as  e = x/T.  If 
the  principles  of  Newtonian  mechanics  are  applied  to  the  falling  body 
after  time  of  impact  where  force  equals  mass  times  acceleration,  then 
the  equation  of  motion  for  the  body  can  be  expressed  as 

Mx*  = -a  A (III -3) 

s c 

Upon  substitution  of  Equation  (III-2)  into  Equation  (III-3) 

A K A Tj 

M2  = - x - “ x (III-4) 

where  x is  measured  relative  to  the  initial  location  of  the  upper  sur- 
face of  the  cushion  and  is  considered  positive  downward.  Rearranging 
the  equation  of  motion  gives 
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A n 


A K 


x + ~ x + -rjz-  x * 0 
MT  MI 

Equation  (III-5)  can  be  cast  in  the  following  canonical  form: 


x + 20  u)  x+cj  x 
c n n 


provided  0^  and  id^  are  defined  as  follows: 


0 = — ns_ 

pc  2a  T (d 


s 


2MT  id 
n n 


- ZS I-  /ST” 

> o 1 ~ v MI 


(II  -5) 


(III-6) 


(HI-7) 


The  parameters  0c  and  (d^  have  physical  significance  as  being  the 
damping  factor  and  natural  frequency  of  the  undamped  system,  respec- 
tively. For  the  case  of  "underdamped"  motion,  the  solution  of  Equation 
(III-6)  can  be  expressed  as 

-0  td  t 

r P n 


x(t)  = e 


c n 


[B^  cos  nt  + B2  sin  ^t] 


where 


•\fl 


(III-8) 


(III-9) 


and  and  B2  are  arbitrary  constants  which  can  be  determined  from  the 
initial  conditions.  Applying  the  initial  conditions 

x(0)  - 0 , 

x;o>  = fig iT  (in-10) 

to  evaluate  the  two  constants  B.  and  B,  gives 

1 L. 

Bx  = 0 


B2  = 


u - p2 
nv  Kc 


(III-ll) 
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Therefore,  the  acceleration  history  can  be  expressed  as 


x(t)  = 


Jn  “n  fc 


(26^  - 1)  sin  at  - 2Q£ 


fl  - ? 


(III-12) 


In  this  research,  only  the  maximum  value  of  the  acceleration  is  of 
interest.  Consequently,  it  is  necessary  to  separate  the  solution  into 
oscillatory  and  decaying  components.  To  accomplish  this,  a trigono- 
metric identity  must  be  introduced  into  Equation  (III-12)  to  obtain  the 
following  expression: 


x(t)  = - 


-Bc  an  t 


1 - 0 


(III-13) 


where 


y = tan 


20  - 1 

c 


^ J 1 - 


(III-14) 


Equation  (III-13)  expresses  the  acceleration  of  the  falling  mass 
as  a function  of  time.  It  is  observed  that  the  character  of  the 
response  is  that  of  a damped  sinusoid  with  the  exponential  and  trigono- 
metric components  of  the  solution  governing  the  transient  behavior. 
Since  the  magnitude  of  the  exponential  or  the  trigonometric  term  can 
achieve  a maximum  value  of  one,  the  coefficient  of  these  terms  governs 
the  absolute  magnitude  of  the  acceleration.  Thus,  in  regard  to  the 
form  of  the  equation,  it  can  be  seen  that 


un/2gh 


1 - 0 


(III-15) 
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Recalling  the  expansion 


2 12  34 

(1-Z)  -1+sZ  +tZ  + . . , , (III-16) 


the  p^ak  acceleration  can  be  expressed  as 


5c  wo)  J 2gh 
max  n v 6 


(III-17) 


where  only  terms  up  through  the  second  order  in  0^  have  been  retained. 
Returning  to  the  primary  variables  o , T,  h,  andQ  , and  Equation  (III-7)  gives 


x ~ g 
max 


/“L\1/2  f1  + I _afs_ 

V“,TJ  [ 8 “SIKJ 


(III-18) 


Inserting  arbitrary  constants  instead  of  the  specific  numerical  coef- 
ficients gives 


.1/2  yl/ 2 , 

V — p 11  --  - 4_  p 11  - - _ 2 

Xmax  1/2  1/2  + U2  3/2  3/2  1 

as  T °s  T 


(III-19) 


where  rj  is  a material  property  dependent  on  temperature. 

Assuming  thermorheologically  simple  behavior  for  the  viscoelastic 


material,  q can  be  expressed  as 


*1(0)  = '10  a(y) 


(III-20) 


where  a(0)  is  the  "shift  factor".  Furthermore,  a(0)  has  been  shown  by 


experience  [ 16]  to  be  expressible  as 


» ”^2  (0”0f))  ” ~C~  0 

a p)  = C1  e U = Cl  e 


(III-21) 


Substituting  into  Equation  (III-20)  gives 


|(9)  - "|Q  e 


it  -C2  u 


(III-22) 


2 "C4 0 

1(0)  = C,  e 4 


(III-23) 
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which  can  be  expanded  In  a Taylor  series  to  give 

ti2(6)  - C3  [i  - C4e  + C5e2]  (in-24) 

where  only  terms  up  through  the  second  order  have  been  retained.  Upon 
substitution  of  this  expression  into  Equation  (111-19),  we  get 

hl/2  hl/2  _ hl/2 

*max  * K1  1/2  _l/2  + K2  “372^372  + K3  3/2  3/2  9 

a t a t a t 

ass 


.1/2 

+ V “ 

\ 3/2  3/2 

a T 


(111-25) 


as  an  expression  relating  the  variables  a , T,  h,  and  6 to  the  peak 

8 

acceleration. 

If  additional  terms  are  retained,  a more  general,  expression  will 
result.  Such  a general  expression  is 


where  the  C ^ (i  * 0,  1,  2,  ...,  N)  and  (j  =■  1,  2,  . . . , M)  are  con- 
stants to  be  determined  by  curve  fitting  procedures. 

Expressions  of  the  form  given  In  Equations  (III-25)  and  (III-26) 


appear  likely  candidates  for  empirically  curve-fitting  cushioning  system 
experimental  data.  This  concludes  Cost's  derivation. 


Verification 

Once  a candidate  relationship  was  selected,  it  was  verified.  A 
statistical  comparison  of  the  G-level  predicted  by  the  function  to  the 
data  base  listed  in  Appendix  A was  used  to  verify  the  model. 

This  was  accomplished  by  modifying  the  Stepwise  Regression  Proce- 
dure given  by  Draper  and  Smith  [ 18] . In  the  Draper  and  Smith  version. 
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variables  are  entered  and  removed  from  the  regression  at  each  stage  on 
the  basis  of  F-tests  performed  on  the  coefficients  of  the  independent 
variables.  This  procedure  was  modified  in  this  research  in  that  the 
variables  were  entered  into  the  regression  equation  and  remain  there  as 
long  as  they  make  a contribution  to  the  overall  correlation.  This  modi- 
fication was  considered  most  appropriate  during  the  preliminary  phases 
of  the  modeling.  It  precludes  the  removal  of  a variable  from  the 
regression  equation  at  one  stage  because  other  variables  can  explain 
most  of  the  variation,  and  then  find  in  a later  stage  that  the  variable 
was  needed  but  was  not  available.  In  the  Draper  and  Smith  version,  the 
establishment  of  critical  F values  can  inadvertently  cause  the  removal 
of  a desirable  variable. 

It  is  considered  most  important  in  the  early  stages  of  model 
development  that  all  variables  that  can  possibly  be  retained  remain  in 
solution.  Once  they  are  discarded  as  irrelevant,  if  is  difficult  to 
reincorporate  them  and  there  is  a substantial  risk  that  they  will  be 
lost.  However,  the  modification  can  cause  the  stages  in  the  solution 
to  carry  superfluous  variables  through  several  iterations  and  ultimately 
retain  them  in  the  final  solution.  This  does  require  some  small  amount 
of  additional  computer  time.  However,  these  unnecessary  variables  are 
easily  identified  and  can  be  discarded  during  the  fine  tuning  of  the 
final  solution  with  little  risk  to  the  validity  of  the  model.  If 
variables  are  discarded  prematurely,  which  could  happen  in  the  Draper 
and  Smith  version  of  the  regression  analysis,  the  validity  of  the  model 
may  be  compromised. 

A program  listing  of  the  Stepwise  Regression  Procedure  that  was 
utilized,  calledMLRD,  is  given  in  AppendixB.  One  form  of  output  from  the 


program  is  printer  plots  in  the  form  of  dynamic  cushioning  curves. 

These  are  a nesting  of  various  thicknesses  of  a particular  cushioning 
at  a particular  drop  height  and  temperature.  The  curves  are  plotted 
using  numbers  corresponding  to  the  thickness  of  cushion.  The  dynamic 
cushioning  curves  that  were  selected  from  the  UAH  study  are  displayed 
in  Appendix  C using  the  MLRD  format.  The  initial  validation  exercise 
for  the  modeling  effort,  then,  consisted  primarily  of  displaying  the 
curve  shape  and  nesting  of  the  model  being  validated  to  determine  if 
the  model  generated  dynamic  cushioning  curves  that  compared  favorably 
with  those  in  the  UAH  study. 

Theoretical  Validity 

The  required  form  of  the  model  as  given  in  Equation  (III-l)  can  be 

identified  in  viscoelastic  theory  as  a phenomenological  constitutive 

equation.  The  model  proposed  by  Cost  was  selected  as  the  most  viable 

candidate,  where  x is  the  maximum  G-level  experienced  by  a body  dur- 
max 

ing  impact  into  a cushion,  and  Cq,  C^,  , Kj*  etc.  are  regression 

constants.  This  equation  was  used  as  the  basic  underlying  structure  of 
the  model  because  it  provides  a functional  relationship  of  all  the 
variables  required  in  Equation  (III-l)  and  has  a theoretical  basis  in 
viscoelastic  theory.  There  was  one  transformation  of  variables  made 
prior  to  using  the  equation.  It  can  be  reasoned  that  when  i)  is  a signif- 
icant aspect  of  the  model,  then  for  0 ** , where  an  exponent  j is  added, 
the  polarity  of  will  oscillate.  To  avoid  this  difficulty,  v was 
transformed  to  °R,  so  the  extreme  temperatures  encountered  in  the  data 


are  always  positive. 
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Validity  of  Initial  Trials 

A program  code  of  Equation  (III-26)  was  run  on  MLRD  and  the  MLRD 
plot  (Figure  8)  shows  that  x was  a progressively  decreasing  value  as 
ag  increased.  This  was  anticipated  to  an  extent  by  Cost  during  his 
finite  element  analysis  [17].  He  generated  a plot  of  peak  acceleration 
versus  static  stress  in  the  region  where  stress  is  proportional  to 
strain  (Figure  9) . 

This  deficiency  in  a aonstitutive  equation  is  not  at  all 
uncommon  as  Meinecke  [15]  suggests  that  even  though  the  behavior  of 
actual  materials  is  usually  different  from  that  predicted  by  various 
classical  theories,  for  engineering  purposes  it  may  be  worthwhile  to 
approximate  the  acvcal  behavior  to  the  idealized  behavior.  But  for 
cases  where  it  is  not  possible  to  approximate  the  actual  behavior  of 
the  material  to  the  idealized  behavior  without  sacrificing  the  accuracy 
in  prediction,  it  is  very  essential  to  consider  the  anomalies  from  the 
idealized  character.  The  departure  of  actual  material  from  its  ideal- 
ized character  may  be  due  to  various  reasons.  For  instance,  the 
deviation  from  the  idealized  character  may  be  such  that  the  stress, 
instead  of  being  directly  proportional  to  the  strain,  is  related  to  it 
in  a complicated  manner.  Solids  display  this  behavior  when  they  are 
stressed  beyond  the  so  called  elastic  range  or  when  the  deformation 
becomes  so  large  as  to  introduce  nonlinearity  between  stress  and  dis- 
placements. Likewise,  liquids  for  which  the  stress  is  nonlinearly 
related  to  rate  of  strain  are  termed  non-Newtonian  liquids.  Equation 
(III-26)  is  a good  representation  of  the  linearity  of  impact  response  but 
a nonlinearity  must  be  introduced  as  an  initial  correction  for  shape  to 
get  a U shaped  curve. 
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Initial  Correction  for  Shape 

All  known  dynamic  cushioning  curves  exhibit  a characteristic  U 
shape.  Consider  a compression  spring  approximately  6 inches  in  diam- 
eter and  a free  height  of  6 inches,  similar  to  those  used  in  the  front 
suspension  of  an  automobile.  On  this  spring  a rigid  steel  plate  that 
serves  as  an  impacting  surface  is  placed.  Then,  a safety  pin,  an 
automobile,  and  a locomotive  are  dropped,  in  turn,  onto  the  spring  from 
approximately  6 inches. 

Now  the  results  of  the  drop  test  are  examined.  The  compression 
spring  is  much  too  stiff  for  the  small  weight  of  the  safety  pin,  so  the 
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spring  does  not  deflect  to  store  the  energy  of  the  fall.  This  rapid 
deceleration  of  the  pin  causes  it  to  experience  high  G-levels.  Because 
the  compression  spring  was  taken  from  an  automobile,  the  spring  deflects 
with  the  fall  of  the  automobile.  Energy  is  stored  by  the  spring  and 
the  load  on  the  automobile  is  proportional  to  the  stiffness  of  the 
spring  and  its  deflection.  For  the  locomotive,  due  to  the  extremely 
heavy  mass,  the  spring  deflects  until  the  coils  of  the  springs  have 
compressed  completely.  Up  to  this  point  the  deceleration  has  progres- 
sed slowly  with  no  high  G-levels.  However,  from  this  point  the  locomo- 
tive will  be  decelerated  to  zero  velocity  almost  immediately  and  will 
experience  high  G-levels. 

Thus,  for  a given  stiffness  and  free  height  of  the  spring,  the 
appropriate  weight  dropped  from  a given  height  will  result  in  the  maxi- 
mum energy  storage  in  the  spring  and  a minimum  G-level  experienced  by 
the  weight.  The  U shapes  of  dynamic  cushioning  curves,  such  as  those 
in  Figure  l,are  the  result  of  this  property  of  bottoming  of  bulk 
cushions.  The  optimal  conditions  exist  at  the  ogive  of  the  curve  as 
indicated  in  Figure  10.  At  this  point  the  maximum  energy  is  stored  in 
the  cushion  with  the  accompanying  minimum  G level. 

A simplified  model  of  a foam  material  such  as  the  one  proposed  by 
Gent  and  Thomas  [ 19]  (Figure  11)  is  helpful  in  visualizing  the  transi- 
tional aspects  of  why  bottoming  occurs.  The  foam  consist  of  thin 
threads  joined  together  to  form  a cubical  lattice.  Mechlin  [20] 
explains  the  bottoming  effect  as  the  result  of  the  ligaments  of  the 
cushioning  material  structure  packing  together,  one  against  another. 


Figure  10.  Typical  dynamic  cushioning  curve  (Minicel,  30-in.  drop  height,  70°F,  1 in.  thickness). 


Figure  11.  Structural  model  of  a foam  material  [22], 


This  situation  could  be  characterized  as  representing  an  entirely  new 
material  with  stress-strain  properties  substantially  different  from  the 
original  cushioning  material.  Accordingly,  the  initial  model.  Equation 
(III-26)  , does  not  address  the  non-linearity  introduced  beyond  the  bot- 
tomed region  of  the  curve,  and  in  itself  is  not  a satisfactory  model  of 
cushion  response.  However,  it  is  reasonable  to  assume  that  the  func- 
tional arrangement  of  the  variables  may  not  be  rearranged  through  the 
bottomed  region  of  the  curve.  This  possibility  was  explored  using  a 
modular  modeling  technique  and  it  was  found  that  the  basic  arrangement 
of  the  variables  has  merit. 


Chapter  IV 


THE  IMPACT  RESPONSE  MODELS 
Development  of  the  General  Model 

The  Initial  model.  Equation  (III-26) , proved  deficient  in  represent- 
ing the  nonlinear  characteristics  of  cushion  response.  However  an 
extensive  literature  search  showed  it  to  be  the  only  known  model  that 
provides  a direct  relationship  of  the  required  variables.  Consequently, 
a modular  technique  suggested  by  Shannon  [ 21]  was  used  to  modify  this 
model  and  construct  a valid  model  of  impact  response.  The  relationship 

of  each  independent  variable  (o  , T,  9 , h)  and  its  effect  upon  the 

8 

dependent  variable  (G-level)  was  studied  and  the  finalized  relationship 
for  each  independent  variable  was  then  incorporated  into  the  model. 

Variable  1.  Drop  Height  (h) 

The  effect  of  drop  height  upon  G-level  that  is  given  in  Equation 
(II-l)  is  based  upon  the  relationship, 

V = /2ih  . (IV-1) 

V is  the  velocity  at  impact  and  is  related  as  the  square  root  of  drop 
height  (h).  Mindlin  [3],  Janssen  [9],  and  others  utilized  this  relation- 
ship of  G-level  versus  drop  height,  and  the  derivation  by  Cost  in 

Equation  (III-26)  is  on  this  same  basis.  It  was  determined  that  drop 

1/2 

height  should  be  incorporated  into  the  model  as  h 
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Variable  2,  Static  Stress  (0  ) 

s 

In  the  UAH  study  [12]  , many  relationships  of  G's  versus  static 
stress  were  investigated  and  it  was  found  that  the  best  agreement  was 
obtained  in  a second  order  polynomial  of  the  natural  log  of  stress. 
Several  polynomials  of  various  orders  of  static  stress  were  tested  in 
this  research  and  a similar  conclusion  was  reached,  namely,  that  a 
second  order  polynomial  was  superior  and  that  the  desired  U shaped 
dynamic  cushioning  curves  of  G’s  versus  static  stress  would  result. 

The  initial  MLRD  Computer  runs  were  made  using  the  following 
relationship: 

G = F(os,  h)  . (IV-2) 

1/2 

The  variables  were  input  with  h and  a second  order  polynomial 

of  0 . The  best  fit  was  obtained  using  the  following  functional 
s 

relationship: 

G - CQ  + C1  h1/2  + C2  h1/2  £n  ag  + C3  h1/2  (in  og)2  . (IV-3) 

These  polynomials  generated  U shaped  curves  similar  to  those  in  the  UAH 
study  (compare  Figure  C-2  and  Figure  12). 

Variable  3,  Thickness  of  Cushion  (T) 

Janssen  [9],  Mindlin  [3],  and  others  suggested  that  G-level  was 
an  inverse  relationship  to  thickness  as  given  in  Equation  (II-2) . This 
seems  intuitively  correct  and  many  attempts  were  made  using  this  rela- 
tionship. T was  introduced  into  these  as  a negative  exponential  such 
-1/2  -3/2  -5/2  -7/2 

as  T , T , T , T , etc.  Computer  runs  showed  good  correla- 
tion with  thickness  input  as  a negative  exponential  of  the  type  given. 
Figure  12  shows  the  typical  nesting  effect  obtained. 


Variable  4.  Temperature  (6) 


Temperature  effects  were  expected  to  be  the  most  difficult  to 
incorporate  and  this  turned  out  to  be  the  case.  It  can  be  reasoned 
that  the  phase  shift  effect  of  temperature  discussed  in  viscoelastic 
theory  for  thermorheologically  simple  materials  produces  the  multi- 
plier effect  of  the  exponentials  shown  by  Cost.  Models  were  tried 
with  several  orders  of  temperature  and  the  polynomial  with  0^,  where 
j * 1,  2,  3,  . . . , n were  the  most  satisfactory.  Figure  13  is  a typical 
plot  using  0^  (where  j = 1,  2,  3)  that  shows  the  shifting  obtained  from 
a cold  temperature  of  -65°F  (C)  through  ambient  (A)  to  hot  160°F  (H) . 

The  General  Model 

Many  relationships  were  tried  and  rejected.  Each  time,  the  basic 
underlying  structure  of  the  variables  was  rejustified  and  new  tormula- 
tions  were  attempted.  The  process  was  repeated  until  a valid  General 
Model  of  impact  response  was  developed.  The  General  Model  is  given 
as  follows: 

0 ■ co + 1 »i/2  S ;u7wo  1 ei  s ciju  <*>  • 

£=Q  k=0  j=l  i=0  (IV -4) 

This  General  Model  incorporates  each  of  the  variables  in  the  man- 
ner prescribed  in  the  modular  modeling  effort.  The  curves  produced  by 
this  model  are  all  U shaped  and  can  be  displayed  using  the  MLRD  plot 
routine.  A series  of  plots  from  the  General  Model  for  various  tempera- 
tures and  drop  heights  are  given  in  Appendix  D for  the  18  and  30  inch 
drop  heights  at  -65°,  70°,  and  160°F.  A comparison  of  these  curves 
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Figure  13.  Concluded 
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demonstrates  the  following: 

A)  The  curves  of  G-level  versus  static  stress  are  U shaped, 
consistent  with  cushioning  theory 

B)  G-level  decreases  as  drop  height  decreases,  consistent  with 
Equation  (II-l) 

C)  The  curves  are  shifted  laterally  with  temperature,  consistent 
with  the  thermorehologically  simple  assumption 

D)  G-level  decreases  as  thickness  increases,  consistent  with 
Equation  (II -2) 

E)  Thickness  curves  are  nested  similar  to  Humbert  and  Handlon  [ 11]  , 
Figure  1. 


Precision  of  the  General  Model 

The  precision  of  this  General  Model,  Equation  (IV-4) , has  to  do 
with  how  well  it  can  be  made  to  represent  a particular  type  of  cushion- 
ing material.  This  can  be  determined  through  comparison  with  an 
experimental  data  base  similar  to  that  in  Appendix  A.  Sensitivity 
analysis  can  be  performed  on  the  model  by  adjusting  the  upper  limits  of 
the  summations,  M,  N,  R,  and  S to  determine  if  it  is  possible  to  obtain 
the  desired  precision. 

A sensitivity  analysis  was  performed  using  the  Minicel  data  of 
Appendix  A.  S is  set  at  1,  R is  set  at  1,  N is  set  at  3,  M is  set  at  2, 
and  the  Minicel  Model  takes  the  following  special  form: 

G-co+S  h"/2  IroTwoS  03  I ci3k*  (^°s)1  • 

i=0  k=0  j=l  i=0  (IV-5) 
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This  is  expanded  to  a 36  term  polynomial  as  follows: 

0 ■ C0  + ffi  [C0100  + C1100  tn  °»  + C2lOo(ta  as)2] 

2 

+ JJ2  [C0200  + C1200  in  °s  + C220o(in  %)  ] 

+ ^L/2  [C0300  + C1300  £n  as  + C230o(in  °s)  ] 

jl/2  n I*  2i 

+ t3/2  LC0111  + Cllll  in  °s  + C2lll(in  as)  J 

h1/2  e2  r \21 

+ t3/2  1_C0211  + C1211  in  as  + C2211 ' n a s'  J 

h 1/2  03  f 2 

+ t3/2  |_C 0311  + C1311  °s  + C231l(^n  °s)  J 

jjl/ 2 g r 2i 

+ t1/2  LC0101  + C1101  in  °8  + C210l(in  °s)  J 

h1/2e2  r 21 

+ jl/2  |_C0201  + C1201  Xn  as  + c220l(in  °s)  J 

h1/2  e3  r 2i 

+ t1/2  [_C0301  + C1301  Xn  °s  + C230l(in  °s)  J 

+ ^3?2  [C0110+  C1110  Xn  as  + C211o(Xn  as)2] 

2 

+ t3/2  LC0210  + C1210  Xn  °s  + C221o(^n  °s)  J 

3 

+ t3/2  [C0310  + C1310  £n  as  + C231o(in  as)2]  * 

(IV-5a) 


45 


Minicel  Model 

To  finalize  the  Minicel  Model,  the  MLRD  program  was  used  with 
Equation  (IV-5),  and  successive  analysis  of  variance  tables  were  con- 
structed as  each  additional  variable  was  entered  into  solution.  These 
ANOVA  tables  were  used  to  determine  if  the  incoming  variable  made  a 
signi  ricar contribution  to  the  overall  correlation.  A Duncan  [22]  F 
test  is  utilized  with  this  format;  it  was  found  that  the  25th  variable 
that  entered  did  not  make  a significant  contribution  to  the  G-level 
response  prediction.  The  test  hypothesis  can  be  stated  as  follows: 

Hq  : the  entering  variable  has  no  effect  upon  the  G-level  response 
prediction. 

: the  entering  variable  has  a significant  effect  upon  the 
G-level  response  prediction. 

The  test  is 


MS 


due  to 


calc 


MS 


(IV-6) 


about 


and  the  null  hyoothesis  can  be  rejected  when  F , > F ...  Utiliz- 

ing  an  a level  of  0.05,  the  final  ANOVA' s and  F tests  (Tables  II  and  III) 
show  that  the  24th  variable  makes  a significant  contribution  but  the 
25th  does  not.  One  further  test  is  made  to  verify  that  all  the  regres- 
sion coefficients  in  the  final  Minicel  Model  are  significant.  This  test 
is  conducted  using  a "t"  statistic  as  follows: 


t = 
n 


C 

n 

S 

n 


(IV-7) 


where 


t is  the  t statistic  for  the  nth  term 
n 


! 

i 

i 

i 

} 
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TABLE  II.  ANOVA  TABLE  AND  F TESTS  FOR  ENTERING  THE  24TH 
VARIABLE  INTO  THE  MINICEL  REGRESSION  EQUATION 


Source 

S.S. 

d.f. 

M.S. 

Fcalc 

F.05 

Decision 

Due  to  all  24 
variables 

5,026,531.3 

24 

209,438.8 

1830.76 

1.52 

5ICN 

Due  to  first  23 
variables 

(5,024,189.8) 

(23) 

(218,443) 

1834.9 

1.52 

SIGN 

Due  to  addition  of 
24th  variable  to 
first  23  variables 

(’341.5) 

(1) 

(2341.5) 

20.6 

3.84 

SIGN 

About  regression 
on  all  24  variables 
(Residual) 

172,858.2 

1511 

114.4 

About  regression 
on  the  first  23 
variables 

175,199.7 

1512 

115.9 

Total 

5,199,389.5 

1535 

TABLE  III.  ANOVA  TABLE  AND  F TESTS  FOR  ENTERING  THE  25TH 
VARIABLE  INTO  THE  MINICEL  REGRESSION  EQUATION 


Source 

S.S. 

d.f. 

M.S. 

Fcalc 

F.05 

Decision 

Due  to  all  25 
variables 

5,026,771.0 

25 

201,070.8 

1758.9 

1.52 

SIGN 

Due  to  first  24 
variables 

(5,026,531.3) 

(24) 

209,438.8 

1830.8 

1.52 

SIGN 

Due  to  addition  of 
25th  variable 

(239.70) 

(1) 

239.70 

2.0 

3.84 

NOT  SIGN 

About  regression 
on  all  25  variables 

172,618.5 

1510 

114.39 

2.0 

About  regression 
on  the  first  24 
variables 

172,858.2 

1511 

Total 

5,199,339.5 

1535 
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C is  the  coefficient  of  the  nth  term 
n 

S is  the  standard  error  of  the  nth  term, 
n 

The  test  hypothesis  is  as  follows: 

Hq  : the  nth  coefficient  is  the  same  as  zero 

: the  nth  coefficient  is  significantly  different  from  zero. 

In  conducting  the  test,  the  null  hypothesis  can  be  rejected  when 

t > t ...  With  ana  level  of  0.05,  the  test  of  all  the  coefficients 
n table  ’ 

in  the  Minicel  Model,  as  given  in  Table  IV,  are  found  to  be  signifi- 
cantly different  from  zero. 

The  resultant  Minicel  Model  is  a 25  term  regression  polynomial. 

The  Minicel  Model  has  a 0.983  correlation  coefficient  which  compares 
favorably  with  the  UAH  data  in  Appendix  A. 

The  reliability  of  the  correlation  coefficient  can  be  tested 
statistically  using  a t statistic  defined  as  follows  [ 23] : 


where 

r = the  correlation  coefficient 

n = the  number  of  samples  used  to  derive  the  regression  equation 

t = the  resulting  number  of  standard  errors  of  r in  the  interval 
n 

between  the  computed  r and  r = 0. 


The  test  hypothesis  is  as  follows; 


H 

H 


0 

1 


r = 0 

r > 0 


TABLE  IV.  TEST  OF  THE  SIGNIFICANCE  OF  THE  REGRESSION  COEFFICIENTS  OF  THE  MINICEL  MODEL 
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In  conducting  the  test,  the  null  hypothesis  can  be  rejected  when 

t > t . This  test  of  the  reliability  of  the  correlation  coeffi- 
n table 

cient  was  conducted  on  the  Minicel  Model . The  null  hypothesis  can  be 
rejected  since  tR  ■ 291.2  > fctable*  The  Minicel  Model  can  be  written 
as  follows: 


. 3.54  ecen  s)2 

15.31  0 h1/2  in  S 

'•39  + ,1/2 

t3/2 

3.34  0 h1/2  (in  Sj2 

. 207.49  0 h1/2 

50.35  0 h1/2 

t3/2 

' t1/2 

T 1/2 

1.43  0 h1/2  (in  Sj2 

6.70  02 in  S 54. 

.66  02  h1'2 

t1/2 

Tl/2 

Tl/2 

11.62  0 2 h1/2in  S 

1.30  03  . 2.08  03  in  S 

1/2 

T 

t1/2  t1/2 

0.23  03  (in  S}2  0.40  03  h1/2  . 0.61  0 

i3  h1/2  in  S 

T1/2 

t3/2 

CM 

CO 

H 

0.09  0J  h1/2  (in  SV 

2 3 1/2 

0.87  0 h ' in 

S 233.77  0 

t3/2  t1/2  t3/2 

28.30  fl  in  S 49.68  r2  26.32  0 2 in  S 6.07  02  (in  S)2 
t3/2  + t3/2  + t3/2  t3/2 


6.15  0 3 in  S 1.O703  (inS)2  3.94  (1V 

t3/2  + t3/2  + t1/2 

Op  ^ 

where  0 = 1 - — and  S = static  stress  in  psi  X 100. 

This  model  can  be  used  to  predict  impact  response  for  Minicel 
cushioning  systems.  The  model  is  expected  to  be  95%  reliable  when  used 
within  the  ranges  of  the  independent  variables  which  are  as  follows: 
h = drop  height  from  12  through  30  inches 


. - 4-  -- 
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a * static  stress  range  from  0.03  to  5.0  psi 
8 

0 * temperature  from  -65°  to  160°F 
T ■ thickness  of  cushion  from  1 through  3 inches. 

The  model  will  predict  with  good  accuracy  at  all  levels  of  the 
independent  variables  within  these  ranges.  Also  it  was  found  that  the 
Minicel  Model  does  a reasonable  job  of  predicting  results  beyond  the 
ranges  stipulated  for  the  independent  variables  as  can  be  seen  when 
the  results  of  tests  of  4 and  5 inch  thick  Minicel  samples  are  compared 
in  Chapter  VI. 

Adjustments  in  Precision 

The  special  form  of  the  General  Model  as  given  in  Equation  (IV-5) 
was  used  in  the  validation  of  the  model  for  the  cross-linked  polyeth- 
ylene foam  Minicel  material  and  gave  a correlation  of  0.983.  If  addi- 
tional precision  had  been  required  the  values  of  S,  R,  M,  and  N could 
have  been  increased  which  may  be  necessary  with  other  materials  but 
gave  only  a minimal  improvement  in  precision  here.  It  was  apparent, 
however,  that  increases  in  M,  which  incorporates  ag  with  exponentials 
over  2 is  in  general  not  worthwhile.  Also,  increases  in  N above  3 that 
incorporate  6 at  exponentials  of  6 over  3 are  of  marginal  value  in 
inproved  precision.  The  changes  in  R and  S that  affect  the  exponen- 
tials of  thickness  and  drop  height  should  be  explored  first  if  addi- 
tional precision  is  required  in  a model  of  a particular  material. 

Once  a special  form  of  the  General  Model  is  found  that  represents 
a particular  cushioning  material,  a measure  of  its  validity  can  be 
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demonstrated  utilizing  the  same  statistical  procedures  as  demonstrated 
for  the  Minicel  Model.  Improvements  in  the  precision  of  the  model  will 
be  reflected  in  increases  in  the  "t"  statistic  associated  with  the  test 
of  the  correlation  coefficient.  Values  comparable  to  those  of  the 


Minicel  Model  are  desirable. 


Chapter  v 


VALIDATION 

The  model  building  process  proceeded  through  many  iterations  of 
development  and  verification  and  culminated  in  the  General  Model  of 
impact  response  stated  mathematically  as 

6 - co + 2 h 1/2  2 7I75T)  2 2 ci4u “s'1  • 

f=0  k=0  j=l  i*0  (V-l) 

This  gives  the  basic  underlying  structure  for  a model  of  impact 
response  for  bulk  cushioning  materials.  The  question  of  how  good  a 
model  has  been  developed  is  answered  in  the  validation  process.  Naylor 
and  Finger  [ 24]  suggest  three  stages  of  validation: 

1)  Verification  of  internal  structure 

2)  Empirical  testing 

3)  Verification  as  a predictor. 

Verification  of  Internal  Structure 

The  validating  process  began  when  the  first  model  of  impact 
response  was  developed.  At  that  time  the  ingredients  of  a model  were 
selected  and  their  relationships  postulated  on  the  basis  of  prior  knowl- 
edge and  existing  theory.  The  basis  of  the  General  Model  was  founded 
in  the  theory  of  viscoelasticity,  and  the  individual  parameters 
(9,  ag,  T,  h,  G)  were  arranged  in  the  model  in  a manner  that  i9  consis- 
tent with  theory  and  intuition.  The  General  Model  that  resulted  from 
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the  modeling  effort  incorporates  the  following  characteristics: 

1)  The  dynamic  cushioning  curves  of  G-level  versus  static  stress 
generated  by  the  model  are  U shaped.  This  is  consistent  with  cushion- 
ing theory. 

2)  The  predicted  3-levels  increase  with  increased  drop  height. 

The  higher  drop  heights  incorporate  more  energy  into  the  system  which 
would  increase  the  e lergy  levels  experienced  by  the  falling  body. 

3)  The  predicted  G-levels  increase  with  reduced  cushion  thickness. 
The  G's  experienced  by  a falling  body  is  dependent  on  shock  pulse  dura- 
tion, Equation  (II-l) , and  a thinner  cushion  would  allow  less  shock 
pulse  time  and  an  accompanying  increase  in  G-levels. 

4)  Temperature  effects  induce  lateral  shifts  in  the  dynamic 
cushioning  curves.  This  is  consistent  with  the  phase  shift  function 
concept  of  viscoelastic  theory.  It  is  also  intuitively  consistent  in 
that  reduced  temperatures  would  be  expected  to  stiffen  the  cushioning 
material  and  require  an  increased  stress  level  for  comparable  shock 
attenuation. 

5)  The  dynamic  cushioning  curves  generated  as  output  from  the 
General  Model  form  a series  of  curves  that  are  nested  within  progressive 
values  of  thickness  and  drop  height  for  all  possible  temperature  condi- 
tions and  static  stress  conditions.  It  can  be  concluded  that  the 
internal  structure  of  the  General  Model  is  intuitively  correct  and  the 
output  of  the  General  Model  is  consistent  with  expectation. 

Empirical  Testing 

The  General  Model  of  Equation  ( V - 1 ) is  hypothesized  as  the  model 
of  impact  response  that  is  applicable  as  the  basic  underlying  structure 
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of  a model  for  any  one  of  the  many  cushioning  materials.  An  impact 
response  model  for  any  one  particular  cushioning  material  can  be  con- 
structed by  establishing  the  summation  levels  M,  N,  R,  and  S,  and  the 
values  of  the  regression  coefficients  in  the  General  Model.  This  can 
be  done  through  a testing  program  that  generates  a data  base  similar  to 
Appendix  A.  Then  an  analysis  program  is  required  that  provides  a least 
squares  fit  of  the  data  base.  The  test  program  required  can  be  similar 
to  the  one  conducted  in  the  UAH  study.  Once  the  test  program  is  com- 
pleted and  the  data  base  established,  an  analysis  must  be  performed  to 
develop  the  model.  The  stepwise  regression  analysis  given  in  Appendix 
B or  a similar  analysis  routine  can  be  utilized. 

When  this  procedure  was  followed  in  constructing  the  Minicel  Model, 
Equation  (IV-8) , the  General  Model  was  used  as  the  basic  underlying 
structure  and  the  statistical  tests  performed  in  verifying  the  Minicel 
Model  serve  to  validate  the  General  Model.  The  ANOVA  Table  for  the 
Minicel  Model  (Table  II),  showing  an  Fca^c  of  1830.76  and  a correlation 
coefficient  of  0.983,  is  indicative  of  the  fit  of  the  Minicel  Model 
to  the  UAH  results.  In  addition  to  the  high  correlation  of  the 
model  with  experimental  data,  the  Minicel  Model  also  demonstrates 
the  characteristic  U shaped  dynamic  cushioning  curves  which  we re  one 
of  the  more  important  aspects  of  the  impact  response  model  considered 
in  the  model  development.  The  dynamic  cushioning  curves  produced  from 
the  Minicel  Model  were  plotted  using  the  MLRD  printer  plot  routine  for 
many  conditions  of  thickness,  drop  height,  and  temperature.  The  clas- 
sic U shape  was  evident  in  all  the  plots,  six  of  which  are  given  in 


Appendix  E . 
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an  additional  measure  of  validity  of  the  model  can  be  seen  when 
the  best  fitting  polynomials  in  the  UAH  study  are  compared  to  the 
Minicel  Model.  The  Minicel  Model  plots  in  Appendix  E have  the  indepen- 
dent variables  at  tne  same  levels  as  the  plots  of  the  UAH  curves  in 
Appendix  C. 

Verification  as  a Predictor 

The  final  test  of  validity  of  the  General  Model  is  to  assess  the 
ability  of  the  model  to  predict  impact  response.  It  was  previously 
demonstrated  that  the  Minicel  Model  does  an  excellent  job  of  predicting 
impact  response  when  compared  to  the  actual  data.  However,  it  must  be 
remembered  that  the  Minicel  Model  is  dependent  upon  these  data  when 
formulating  its  predictions.  It  is  a limited  dependency  in  that  it  is 
using  all  2709  data  points  in  predicting  a particular  dynamic  cushion- 
ing curve,  and  in  fact  only  a small  portion  of  the  data  base,  approxi- 
mately 50  data  points,  apply  directly  to  the  particular  conditions 
with  the  independent  variables  at  the  appropriate  levels.  The  UAH  best 
fitting  polynomials,  however,  are  derived  using  only  those  data  points 
where  the  independent  variables  are  at  the  appropriate  levels. 

To  fully  verify  the  model  as  a predictor,  a data  base  independent 
of  that  used  to  generate  the  model  must  be  utilized.  Three  such  data 
bases  are  available:  a data  base  of  1,  2,  and  3 inch  thickness  Minicel 

material  at  -65°,  70°,  and  160°F  and  at  27  inch  drop  height,  and  a data 
base  of  4 and  5 inch  Minicel  material.  The  27  inch  drop  height  data 
are  given  in  Appendix  F and  the  4 and  5 inch  data  in  Appendix  G.  The 
27  inch  data  are  contained  within  the  range  of  the  independent  variables 
and  was  not  used  in  formulating  the  Minicel  Model.  However,  the  4 and 
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5 inch  data  are  beyond  the  data  extremes  of  the  developed  model.  Also 
the  4 and  5 inch  samples  themselves  were  not  homogeneous.  The  samples 
of  the  4 inch  material  were  tv/o-piece  cushions  which  were  2 inches 
thick.  The  5 inch  material  was  made  using  a 2 inch  and  a 3 inch  piece. 
This  stacking  is  representative  of  how  cushioning  is  actually  used  in 
shock  isolation  systems  requiring  thicker  sections  than  the  maximum 
manufactured  thickness  of  3 inches.  Whether  the  Minicel  Model  can  do 
an  adequate  job  of  predicting  the  impact  performance  of  these  stacked 
samples  is  questionable.  The  model's  ability  to  predict  adequately 
under  these  circumstances  would  indi:ate  that  the  stacking  did  not 
significantly  perturb  the  cushioning  performance  from  that  encountered 
in  the  1,2,  and  3 inch  continuous  samples  that  form  the  basis  of  the 
model . 

To  determine  statistically  how  well  the  model  fits  a set  of  inde- 
pendent data.  Box  and  Draper  [ 25]  suggest  it  is  appropriate  to  investi- 
gate the  bias  and  variance  of  the  predictor.  Two  statistical  tests  can 
be  utilized  for  this  purpose.  One  test,  a test  of  means,  determines 
whether  there  is  a bias  in  the  predicted  values  of  the  model  when  com- 
pared to  actual  data  values.  The  other  test,  a test  of  variances, 
determines  if  the  variations  of  the  predicted  values  inherent  in  the 
model  are  comparable  with  the  variations  in  experimental  values. 

Test  of  Means 

In  all  the  data  bases,  including  the  UAH  data,  there  are  three 
replications  of  the  same  conditions  of  the  independent  variables.  These 
three  replications  can  be  considered  a cell.  Then,  the  mean  of  the 
G-levels  of  the  three  data  values  in  a cell  (G^,  i = 1,  2,  3)  can  be 
compared  with  the  predicted  G-level  from  the  model  for  that  cell,  and 
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it  is  reasonable  to  expect  the  sum  of  the  differences  to  vanish.  Any 
difference  that  cannot  be  reasonably  attributed  to  s.impling  error  must 
be  considered  a bias  that  is  introduced  because  the  mod?!  values  are 
not  good  predictors. 

The  first  step  in  this  test  is  to  formulate  a difference  between 

the  data  values  in  a cell  and  the  model  prediction  for  that  cell.  This 

can  be  expressed  mathematically  as  follows: 

a.  = (CM.  - GD.)  (V-2) 

J J J 

where 

j = number  of  a cell  with  fixed  values  of  a,  os>  T,  h 

A.  = cell  difference  for  cell  j 

GMj  = the  G-level  predicted  by  the  model  for  cell  j 

GD.  = the  mean  value  of  G-level  for  the  three  data  values  in 
J 

cell  j 


GD. 

J 


Then  S,  rhe  standard  deviation  of  all  the  cells,  can  be  defined 
as  follows: 


S = 


(V-3) 


where 


N = the  number  of  cells  in  the  data  base 


G = the  grand  mean  of  all  the  cell  differences,  G = 


I 
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The  hypothesis  to  be  tested  is  based  on  the  expected  value  of  the 
differences  which  can  be  written  E(A^)  and  stated  as  follows: 

HQ:  E(A.)  = 0 

Hl:  E(V)  * 0 

A two  tailed  "t"  test  is  used  '..'here  the  t statistic  is  given  as 


_ G /n~ 
m s 

where  t = the  test  statistic  for  the  model. 


(V-4) 


m 

The  test  compares  t^  with  anc*  t^ie  nuH  hypothesis  can  be 

rejected  when  t '>  t ^ . Rejection  of  the  null  hypothesis  implies 
that  E(Aj)  is  sufficiently  greater  than  zero  as  to  be  unexplainable  as 
sampling  error  and  therefore  the  model  predictions  would  appear  not  to 
be  representative  of  the  data. 


Test  of  Variances 

The  other  test  of  goodness  of  fit  of  the  Minicel  Model  with  actual 
data  determines  how  the  variations  in  the  prediction,  using  the  model, 
compare  with  the  variations  in  the  data.  In  each  data  base,  the  data 
points  in  each  cell  are  the  replications  of  the  experiment  for  each  set 
of  conditions  of  stress  level,  drop  height,  and  temperature.  The  sum 
of  these  variations  for  each  cell  can  be  written  as  follows: 


R 

T.  = £ (Gi  " ^j)2  (V"5) 

i=l 


where 

i = the  number  of  the  sample  in  the  cell  being  considered 
T.  = the  sum  of  squares  of  the  variation  for  cell  j 
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R = the  number  of  replications  in  each  cell 

= G-level  value  of  the  data  point  being  considered 
GDj  = mean  value  of  the  data  values  in  the  cell  being 
considered  from  Equation  (V-2) . 

2 

Then  o,  can  be  defined  as  the  data  within-cells  variance  which  is  found 
d 

by  summing  the  variations  within  all  the  cells: 


2 

°d 


N 


izl 

N X df 


(V-6) 


where  df  = degrees  of  freedom  in  each  cell. 

The  hypothesis  to  be  tested  is  whether  the  variance  of  the  data 
base  is  equal  to  the  variance  in  the  model  predictions. 

The  test  hypothesis  can  be  expressed  as  follows: 


H 

H 


0 : 
1 : 


where 

2 

ad  = the  variance  in  the  data  as  given  in  Equation  (V-6) 

2 

= the  variance  in  the  predicted  values  from  the  model. 

An  "F"  value  is  used  to  test  the  ratio  of  variances  and  the  F statistic 
is  defined  as  follows: 


calc 


(V-7) 


The  model  variance  is  computed  during  the  regression  procedure 
as  the  Residual  Mean  Square  in  the  analysis  of  variance  of  the 
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model  being  tested.  The  test  compares  Fca^c  with  ^ca^c 

computed  from  Equation  (V-7)  and  Fta^^e  is  set  using  a = 0.05  and  is 
tabled  according  to  the  degrees  of  freedom  in  the  numerator  (the 
degrees  of  freedom  in  the  data  base  being  considered)  and  the  degrees 
of  freedom  of  the  demoninator  (the  degrees  of  freedom  of  the  model 
being  tested) . 

Rejection  of  the  null  hypothesis  indicates  that  there  is  a dif- 

2 2 

ference  between  the  model  variance,  a,,  and  the  data  variance,  o , that 

* d*  ’ ITT 

cannot  be  attributed  to  sampling  error.  This  implies  that  there 
is  a significant  difference  in  the  variance  of  the  model  when  com- 
pared with  the  actual  data  and  that  the  model  is  not  representative  of 
the  data. 


Prediction  Test  Results 

The  test  of  means  and  variances  were  conducted  on  the  Minicel 
Model,  Equation  (IV-8) , using  the  27-inch  drop  height  data  in  Appendix 
F and  the  4 and  5 inch  data  samples  of  Appendix  G.  The  results  are 
given  in  Table  V and  show  the  following: 

1)  27 -inch  drop  height  data  - The  27 -inch  drop  height  data  are 
within  the  range  of  the  independent  variables  used  in  the  Minicel  Model 
and  the  null  hypothesis  cannot  be  rejected  in  the  test  of  means  or  test 
of  variances.  This  appears  to  indicate  that  there  is  no  significant 
difference  between  the  means  and  variances  of  the  data  and  the  values 
predicted  by  the  model.  The  model  appears  to  be  a statistically  valid 
predictor  of  these  data. 

2)  4 and  5 inch  data  - The  test  results  for  the  4 and  5 inch 
data  show  that  the  Minicel  Model  gives  good  predictions  of  the 
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data  means  at  the  4 inch  thickness.  The  test  of  means  on  the  4 inch 
showed  no  significant  difference  but  the  test  on  the  5 inch  showed  a 
significant  difference  \;hich  is  not  surprising  because  the  model  is 
projecting  2 inches  beyond  its  range  at  the  5 inch  thickness.  The  test 
of  variances  shows  that  the  model  variance  at  the  4 inch  thickness  is 
significantly  larger,  but  this  is  not  true  at  the  5 inch  level.  The 
model  dispersion  is  comparable  to  the  dispersion  of  the  5 inch  data. 

There  were  reservations  as  to  whether  the  model  could  predict  the 
4 and  5 inch  materials  since  the  samples  were  stacked  and  not  homoge- 
nous and  since  4 and  5 inches  are  beyond  the  ranges  of  the  independent 
variables  used  in  the  Minicel  Model,  The.  results  of  the  tests  of  means 
and  variances  show  that  the  model  was  not  completely  acceptable  at 
either  the  4 inch  or  5 inch  thickness.  However,  the  tests  indicate 
that  model  means  were  comparable  with  the  data  at  the  4 inch  level  and 
the  variances  were  comparable  at  the  5 inch  level. 
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Chapter  VI 


OPTIMIZATION 

The  procedure  which  was  designed  in  this  research  to  perform  cushion 
system  optimization  can  best  be  described  as  a constrained  sequential 
search  technique.  The  technique  uses  the  knowledge  that  all  dynamic 
cushioning  curves  are  U shaped,  and  searches  for  limits  of  acceptable 
G-level  values  along  these  curves.  The  procedure  is  a direct  search 
technique  in  that  the  mathematical  model  of  a material  such  as  the 
Minicel  Model,  Equation  (IV-8)  , is  used  directly  as  the  objective  func- 
tion in  the  optimization  procedure. 

The  first  step  in  the  optimization  exercise  is  to  identify  the 
parameters  in  an  optimal  cushioning  system  design;  then,  the  objective 
or  goal  of  the  optimization  procedure  must  be  identified.  The  con- 
straints on  the  procedure  and  the  objective  function  are  formulated 
and  finally  the  outputs  from  the  optimization  routine  itself  are 
identified. 


Formatting  for  Optimal  Cushion  System  Design 
An  optimal  cushioning  system  is  a system  that  provides  the  neces- 
sary shock  isolation  to  the  protected  item  at  a minimum  cost.  Because 
the  cost  of  a cushioning  system  is  dependent  upon  the  amount  of  cushion- 
ing material,  the  optimal  cushioning  system  will  employ  the  minimum 
thickness  of  cushion  needed.  The  optimal  point  on  a dynamic  cushioning 
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curve  such  as  that  shown  in  Figure  10  is  the  minimum  G's  possible  which 
provides  the  maximum  protection  for  that  particular  thickness  of 
cushioning.  Therefore,  it  is  the  identification  of  this  optimal  point 
that  must  be  determined  in  selecting  the  optimal  cushion  system. 

A different  dynamic  cushioning  curve  is  required  to  depict  cushion 
performance  for  each  drop  height,  thickness,  and  temperature  of  the 
cushion;  the  optimal  point  is  different  for  each  curve.  When  one  con- 
siders the  many  types  of  cushioning  materials  and  a different  curve  for 
each  condition,  it  can  be  seen  that  a very  large  library  of  curves  is 
required  to  present  even  the  most  likely  conditions  of  drop  height, 
temperature,  and  thickness  of  a few  candidate  cushioning  materials. 

A model  of  impact  response  for  a cushioning  material,  e.g.  the 
Minicel  Model,  computes  the  impact  response  directly  from  the  values  of 
the  independent  variables.  This  precludes  the  need  for  a library  of 
dynamic  cushioning  curves  in  predicting  impact  response.  It  is  also 
possible,  using  the  cushioning  models  for  various  materials,  such  as 
the  Minicel  Model,  to  determine  an  optimal  cushioning  design  for  each 
material,  if  such  a design  exists.  It  is  also  important  to  provide  the 
designers  with  the  maximum  amount  of  useful  design  information  in  the  output. 

Specification  of  the  Cushion  System  Constraints 

The  development  of  a valid  mathematical  model  of  impact  response 
for  a particular  cushioning  material  provides  a vehicle  for  the  applica- 
tion of  optimization  techniques.  This  model  can  be  used  as  an  objective 
function  and  explored  for  an  optimal  cushioning  system.  This  can  be 
done  by  defining  cushioning  system  design  requirements  in  terms  of  the 
external  environment  and  the  amount  of  exposure  the  protected  item  can 
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The  specification  of  the  external  environment  must  contain  a 

definition  of  those  aspects  of  the  environment  that  have  an  Impact  upon 

cushion  system  design.  This  would  -include  a quantitative  measure  of 

the  magnitude  of  the  maximum  shock  pulse  to  which  the  system  will  be 

exposed,  which  is  usually  specified  in  terms  of  drop  height.  Also 

since  temperature  has  a significant  effect  on  impact  response,  the 

specification  of  the  external  environment  should  include  the  range  of 

temperatures  to  which  the  system  will  be  exposed  (0  , , 0 ) . 

min  max 

The  specification  of  the  survivability  of  the  protected  item  is 

given  in  terms  of  its  ability  to  withstand  shock.  The  maximum  shock 

pulse,  in  G's,  that  the  item  can  withstand  is  given  as  the  fragility 

level  of  the  item  (GL  ) . 

max 

These  considerations  concerning  the  G-levels  and  temperatures  are 
incorporated  into  the  model  as  input  values  that  are  utilized  to  con- 
strain the  optimal  search  to  feasible  solutions. 

Construction  of  the  Objective  Function 

Consider  the  parameters  of  the  general  cushion  design  model  that 
can  be  expressed  mathematically  in  functional  notation  as 

G = F(os,  T,  0,  h)  . (VI-1) 

For  the  purpose  of  optimal  design,  GL  is  identified  as  the  fragility 
level  of  the  protected  item.  The  temperature  parameter,  0,  must  con- 
sider the  range  of  temperatures  of  the  external  environment  0^  to 
^max*  '^ie  suPerimPose<*  dynamic  cushioning  curve  format,  shown  in 
Figure  2,  can  be  used  for  this  purpose  wherein  the  curves  of  extremes 
of  the  temperature  range  are  superimposed  upon  the  ambient  curve. 
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One  parameter  in  the  model,  drop  height  (h) , can  be  determined 

based  on  published  testing  standards.  Other  parameters,  such  as  the 

fragility  level  of  the  item  to  be  protected  (GL^^)  and  the  temperature 

range  (0  . , 0 ) , are  usually  specified  by  overall  system  require- 

min  max 

ments.  For  example,  if  the  cushioning  system  is  to  be  designed 
for  a military  container  for  a missile,  the  fragility  level  of  the 
protected  item  (GLmqv)  will  be  specified  as  part  of  the  hardware  speci- 
fications in  the  missile  system  design  criteria.  The  temperature  range 
^min*  ^max^  t^e  ®i88He  system  requirements  and  the  drop 

height  is  specified  in  various  Military  Standards  depending  on  total 
container  weight.  These  parameters,  GL  , 0 . , 0 , and  h,  are  the 

exogenous  variables  in  the  optimization  model.  Equation  (VI-1)  can 
now  be  written  as  follows: 

a = F[o  ,T,  (0  ,0  ) , h]  (VI-2) 

max  1 s * min’  max7 » 1 ' 7 

where  GL  , (0  . , 0 ),  and  h are  inputs  to  the  equation  which  are 

determined  by  the  cushion  system  design  requirements.  The  optimization 

technique  will  search  out  the  optimal  design  from  this  expression,  and 

the  optimal  design  solution  will  be  expressed  in  terms  of  a and  T. 

s 

The  optimal  design  searches  upon  the  functional  relationship 
expressed  in  Equation  (VI-2)  and  with  the  inputs  introduced,  reduces 
to  the  selection  of  the  minimum  thickness  of  cushion  that  will  perform 
satisfactorily  at  a static  stress  condition  that  is  determined  in  the 
search.  Classical  search  techniques  would  establish  Equation  (VI-2) 
as  the  model  to  use  in  the  search  procedure  and  it  would  be 
rearranged  with  T as  the  dependent  variable  and  the  objective  function 
as  follows: 


min 


- 


GL 


max 


0,  h] 
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The  search  would  be  subject  to  certain  constraints  such  as 

[0  , < 9 < 9 ] and  [GL  < GL  ] , and  a search  made  for  T . using 

min  max  max  min 

one  of  many  search  techniques. 

However,  it  is  not  apparent  how  Equation  (Vl-2)  can  be  rearranged 
mathematically  into  the  simple  form  of  Equation  (VI-3)  where  T is  solved 
for  directly,  nor  is  it  necessarily  desirable  to  do  so.  It  is  important 
to  recognize  that  when  the  exogenous  variables  are  introduced  into  the 
objective  function,  a three-dimensional  search  on  og  and  T is  all  that 

is  required.  Since  the  general  shape  of  the  curve  of  G versus  a is 

s 

known,  i.e.,  it  is  a classical  U shaped  dynamic  cushioning  curve,  Equa- 
tion (VI-2)  can  be  searched  directly  with  little  difficulty  and  with 
reasonable  efficiency. 

Once  Equation  (VI- 1)  is  rewritten  as  Equation  (VI- 2)  with  the 

values  of  h,  (0  . 9 ),  and  GL  introduced,  we  have  an  objective 

* min  max  * max  J 

function  that  can  be  written  as  follows: 

GL  = F ( o , T,  (0  . , 0 ),  h) 

max  \ s min  max  / 

where  GL  ,0,0  , and  h are  input  constraints.  This  objective 

max  min  max  J 

function  is  then  searched  using  the  direct  search  routine. 

The  Direct  Search  Routine 

The  direct  search  routine  (CUSHION  OPT)  is  a three  stage  process. 
The  initial  stage  involves  the  selection  of  a material  type.  It  is 
anticipated  that  a data  base  will  eventually  be  available  that  contains 
valid  impact  response  models  for  many  of  the  different  types  of  bulk 
cushioning  materials.  The  model  of  each  type  of  material  in  the  data 
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base  can  be  investigated  by  the  direct  search  routine  to  determine  the 
feasibility  of  meeting  the  design  objective  function.  The  optimal 
design  conditions  of  minimum  thickness  and  acceptable  static  stress 
range  will  be  output  for  each  type  of  cushion  in  the  data  base  if  an 
optimal  design  exists. 

The  initial  step  is  to  select  the  type  of  material.  Then,  a search 
is  initiated  on  a minimum  thickness  cushion  (1  inch)  at  the  drop  height 
and  minimum  temperature  (Qmjn)  stipulated  in  the  design  requirements. 
These  values  are  input  into  the  objective  function  and  the  search  is 
conducted  across  the  static  stress  spectrum  at  the  temperature  extreme, 
&min  an<*  then  9mavt  to  determine  the  feasibility  of  meeting  the  stipu- 
lated GLmax>  It  is  also  necessary  that  the  acceptable  static  stress 
range  at  GLmax  be  greater  than  0.2  psi  to  permit  design  flexibility 
and  preclude  creep  problems  within  the  design.  (This  0.2  psi  value  can 
be  changed  at  the  discretion  of  the  designer.). 

If  the  search  of  the  acceptable  static  stress  range  on  the  minimum 
cushion  thickness  is  successful,  the  answer  is  output  as  a feasible 
solution.  If  it  is  not  successful,  the  thickness  is  incremented 
1/2  inch  and  the  search  repeated.  A flow  chart  of  the  CUSHION  OPT 
search  routine  is  given  in  Figure  14  and  the  program  code  is  given  in 
Appendix  H . 


Results 

The  CUSHION  OPT  Program  output  is  given  in  .i  e form  of  superimposed 
dynamic  cushioning  curves  such  as  the  typical  one  given  in  Figure  15 . 

In  Figure  15,  the  design  objective  function  had  the  temperature  range 
of  -65°  to  160°F  and  a 30-inch  drop  height.  The  fragility  level  was 
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Figure  14  . Flow  chart  of  optimization  direct  search  routine  (CUSHION OPT) 
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50  G's  which  is  indicated  with  a dotted  horizontal  line.  The  optimal 
thickness  is  determined  to  be  2.5  inches  and  the  feasible  static  stress 
range  is  found  to  be  0.60  to  1.20  psi.  One  of  the  inherent  advantages 

of  using  math  models  is  seen  here  in  the  capability  of  determining 

G-level  response  at  thicknesses  other  than  those  in  the  data  base.  In 
this  instance  a 2.5-inch  thickness  of  cushion  was  optimal;  tests  were 
not  conducted  at  this  thickness  in  the  data  base.  This  advantage  can 
also  be  seen  in  Figure  16  which  is  another  output  of  the  optimization 
program  where  the  temperature  range  has  the  non-data  base  values  of 
-20°  and  120°F.  The  optimal  thickness  is  2.0  inches  and  the  stress 

range  is  0.55  to  0.85  psi  to  give  50  G protection.  Comparison  of 

Figure  15  with  Figure  16  demonstrates  the  effect  of  temperature  on 
optimal  cushion  design.  When  the  temperature  range  was  relaxed  from 
the  extremes  of  -65°  through  160°F  to  -20°  through  120°F  and  all  other 
exogenous  variables  kept  the  same,  the  thickness  of  cushion  required 
for  50  G protection  dropped  from  2.5  inches  to  2 inches.  Another  com- 
parison can  be  made  between  Figure  15  and  Figure  17,  which  demonstrates 
the  increase  in  thickness  of  cushion  required  when  the  fragility  level 
of  the  protected  item  is  lowered  from  50  G's  to  40  G's.  One  other 
comparison  can  be  made  between  Figure  15  and  Figure  18,  which  demonstrates 
how  a reduced  drop  height  requirement  reduces  the  thickness  of  cushion 
required. 

One  additional  advantage  in  using  the  superimposed  dynamic  cushion- 
ing curve  format  for  the  output  format  is  that  the  designer  is  presented 
with  a convenient  tool  for  minimizing  the  creep  tendency  of  the  mate- 
rials by  maintaining  as  low  a stress  condition  as  possible  and  yet  be 
able  to  ascertain  the  response  of  the  stress  level  selected. 
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Figure  17.  Optimal  design  output,  3 in.  Minicel. 
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Chapter  VII 


CONCLUSIONS  AND  REOOM4ENDATIONS 

There  has  been  limited  use  made  of  bulk  cushioning  in  shock  isola- 
tion systems  for  commercial  and  military  equipiaent.  One  reservation 
that  designers  have  had  regarding  bulk  cushioning  is  an  inability  to 
predict  cushioning  performance  at  temperature  extremes.  Consequently, 
many  designs  incorporate  mechanical  suspension  systems  such  as  elasto- 
meric mounts  or  springs  and  dash  pots  and  other  types  of  shock  isola- 
tion systems  that  are  bulky,  very  heavy,  require  substantial  structural 
interfaces,  and  are  very  expensive. 

Oue  result  of  this  research  is  the  development  of  a valid  model  of 
bulk  cushioning  response  that  provides  a basis  for  improving  the  pre- 
dictability of  temperature  performance  for  bulk  cushioning  materials. 
Resultant  increases  in  the  use  of  bulk  cushioning  for  shock  mitigation 
systems  will  generate  dollar  savings  in  the  accompanying  reduced  pro- 
curement and  development  costs. 

Conclusions 

The  developed  General  Model  of  the  impact  response  of  bulk  cushion- 
ing materials  is 

C"  C0  + I h#2  I ^uTSk)  I eJ  I • 

i=0  k=0  j=l  i=l  (VII-1) 
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The  model  Is  predicated  on  viscoelastic  theory  and  incorporates 
the  effect  of  temperature,  stress,  drop  height,  and  thickness  of 
cushion  upon  the  impact  response  of  a cushioning  system.  This  General 
Model  provides  the  basic  underlying  structure  of  impact  response  of 
any  one  of  the  many  bulk  cushioning  materials  used  for  shock  isolation. 

A sensitivity  analysis  can  be  run  on  the  values  of  S,  R,  N,  and  M to 
obtain  the  precision  desired  for  an  impact  model  of  a particular 
cushioning  material. 

Models  that  are  predicated  on  the  basic  underlying  structure  of 
the  General  Model  are  better  predictors  of  impact  response  than  the 
dynamic  cushioning  curves  currently  being  utilized,  because  the  effect  of 

temperature  has  been  incorporated  into  the  model.  The  Minicel  Model  is 

3 

one  such  model  that  was  constructed  for  the  Hercules,  Inc.  2 lb/ft 
Minicel  material  using  the  General  Model  as  the  basic  underlying 
structure.  The  Minicel  Model  is  a 25-term  polynomial  given  in  Equa- 
tion (IV-8) . The  correlation  of  the  Minicel  Model  with  actual  data 
demonstrates  the  validity  of  the  models  and  their  value  as  predictors 
of  impact  response.  The  Minicel  Model  showed  high  correlation  with 
the  actual  data  within  the  ranges  of  the  variables  and  also  showed 
promise  as  a predictor  beyond  the  variable  ranges. 

The  development  of  a valid  model  of  impact  response  of  bulk 
cushioning  materials  defines  the  relationships  and  interrelationships 
of  the  variables  in  response  to  impact.  Using  this  basis,  it  war  pos- 
sible to  employ  a search  technique  (CUSHION  OPT)  to  determine  optimal 
cushion  design  and  display  the  findings  in  terms  of  superimposed  dynamic 
cushioning  curves.  The  functional  form  of  the  model  provides  the 
advantage  of  determining  impact  response  at  non-tested  levels  of  the 
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variables  with  confidence  and  precludes  the  need  for  a library  of 
dynamic  cushioning  curves  for  all  the  different  combinations  of 
conditions. 

Once  a valid  model  of  a particular  cushioning  material  has  been 
developed,  it  can  be  incorporated  into  the  CUSHION  OPT  optimization 
program.  This  program  accepts  the  design  requirements  for  a shock 
isolation  system  and  computes  and  provides,  in  the  form  of  superimposed 
dynamic  cushioning  curves,  the  optimal  design  for  each  cushioning 
material  in  the  data  base,  if  one  exists.  The  outputted  superimposed 
dynamic  cushioning  curves  give  the  pertinent  information  needed  in  the 
optimal  design  of  a shock  isolation  system. 

Recommendations 

It  is  recommended  that  the  model  of  impact  response  of  bulk 
cushioning  materials.  Equation  (VII-1)  , be  used  as  the  basic  under- 
lying structural  design  for  cushioning  systems.  This  model  is  consider- 
ably better  than  any  previous  basis  of  design.  The  optimization  pro- 
gram used  in  conjunction  with  the  model  can  be  used  to  provide  accurate 
predictions  of  shock  mitigation  system  performance  in  a time  saving 
manner  and  in  a useful  format.  It  is  rea.  -mable  to  expect  that  in  the 
design  of  shock  isolation  systems,  considerable  savings  can  be  realized 
in  design  time  and  cost  savings  by  using  these  more  accurate  response 
predictions  prior  to  prototype  fabrication  and  test. 

Further,  it  is  recommended  that  the  optimal  design  program,  CUSHION 
OPT,  be  used  to  generate  the  optimal  design  of  bulk  cushioning  systems. 
Cushioning  system  designers  that  have  a large  scale  digital  computer 
capability  should  be  encouraged  to  utilize  this  procedure.  As  models 
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of  additional  cushioning  materials  become  available,  CUSHION  OPT  is 
provisioned  to  permit  updating  to  incorporate  additional  design  alter- 
natives using  the  new  materials.  To  insure  maximum  utilization  of  the 
improved  predictive  capability  of  the  General  Model,  it  is  further 
recommended  that  superimposed  dynamic  cushioning  curves  of  the  most 
likely  conditions  of  the  independent  variables  ( , og,  h,  T)  be  pub- 
lished and  made  available  to  cushioning  system  designers  who  do  not 
have  access  to  a computer  facility.  It  is  recognized  that  one  of  the 
major  advantages  inherent  in  the  CUSHION  OPT  program,  that  of  obtaining 
the  exact  design  constraints  needed,  becomes  inoperative  and  therefore, 
manual  interpolation  will  be  necessary  when  using  published  curves. 
However,  the  difficulties  associated  with  using  and  maintaining  a 
library  of  superimposed  dynamic  cushioning  curves  are  warranted  when  the 
savings  accrued  through  the  use  of  the  improved  predictions  are 
considered. 

Additional  drop  test  programs , similar  to  the  one  conducted  on 
Minicel,  should  be  conducted  on  the  bulk  cushioning  materials.  A model 
of  impact  response  similar  to  the  Minicel  Model,  Equation  (IV-8),  using 
the  General  Model  as  the  basic  underlying  structure  should  be  con- 
structed for  each  new  material.  This  model  can  then  be  entered  into 
the  CUSHION  OPT  program  to  provide  an  additional  optimal  design  alter- 
native in  terms  of  this  new  material.  It  is  further  recommended  that 
the  drop  test  programs  on  additional  materials  be  conducted  using  levels 
of  the  independent  variables  (t  , o , h,  T)  that  bracket  the  values  of 
the  variables  that  might  occur  in  cushion  system  design.  This  increase 
in  the  range  of  the  variables  would  insure  that  the  model  is  predicting 
values  within  the  range  of  the  model.  Particular  attention  should  be 
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given  to  insuring  that  the  range  of  thickness  is  sufficient  to  obtain 
G-levels  as  low  as  10  G’ s,  which  is  not  uncommon  in  some  of  the  more 
fragile  optical  and  electronic  hardware.  Also,  since  only  the  lower 
portions  of  dynamic  cushioning  curves  are  used  in  optimising  cushion 
design,  the  test  programs  can  be  abbreviated  in  the  understressed  and 
over stressed  regions. 

It  is  also  reconmended  that  the  General  Model  of  impact  response 
be  considered  as  a basis  for  future  research.  The  mathematical  formu- 
lation of  impact  response  should  provide  a vehicle  to  be  utilized  for 
the  rigorous  analysis  of  impact  response.  For  example,  one  particularly 
lucrative  area  is  that  the  model  be  used  as  a phenomenological  constitu- 
tive equation  in  advancing  the  viscoelastic  theory  of  bulk  cushioning 
materials. 
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Appendix  A 

UAH  DATA,  12,  18,  24  AND  30  INCH 


STRESS  LEVELS  (PSI) 


CORRELATION  COEFFICIENTS  OF  THE  BEST  FITTING 
POLYNOMIALS  IN  THE  UAH  STUDY 


Temperature 

(°F) 

Drop  Height 

(in.) 

12  in. 

18  in. 

24  in. 

30  in. 

1 

0.97 

0.98 

0.98 

0.97 

-65 

2 

0.99 

0.99 

0.96 

0.99 

3 

0.99 

0.98 

0.99 

0.98 

1 

0.98 

0.98 

0.99 

0.98 

70 

2 

0.98 

0.96 

0.97 

0.98 

3 

0.97 

0.98 

0.97 

0.97 

1 

0.98 

0.98 

0.97 

0.98 

160 

2 

0.97 

0.97 

0.94 

0.97 

3 

0.93 

0.98 

0.96 

0.97 

NOTE:  This  table  gives  the  sample  correlation  coefficients 

from  an  analysis  of  regression  variance  between  the 
data  and  a polynomial  fit  of  the  form 

2 

y.  = bn  + b^nx.  + b„(£nx.)  . 

1 0 1 i 2 i 


SELECTED  BEST  FITTING  POLYNOMIALS  IN  THE  UAH  STUDY 
(Hercules  Minicel,  2 lb/ ft2  Density) 


Thickness 

Temperature 

(in.) 

(°F) 

Design  Curve  Equation 

12  in.  Drop  Height 

-6; 

y = 377.74  - 142.48  in  x + 14.78  (in  x)2 

1 

70 

y = 278.24  - 118.34  in  x + 14.44  (in  x)2 

160 

y = 197.11  - 84.21  in  x + 11.27  (in  x)2 

-65 

y = 367.89  - 131.51  m x+  12.22  (in  x) 2 

2 

70 

y = 201.97  - 78.34  xn  x + 8.34  (in  x)2 

160 

y = 142.03  - 55.43  in  x + 6.31  (in  x)2 

-65 

- — j 

y = 329.67  - 118.13  in  x + 10.86  (in  x)“ 

3 

70 

y = 159.93  - 58.41  in  x + 5.77  (i.n  x)  “ 

160 

y = 105.43  - 36.66  in  x + 3.73  (in  x)  “ 

24  in.  Drop  Height 

-65 

y = 691.02  - 301.76  in  x + 36.16  (in  x)2 

1 

70 

y - 403.76  - 193.48  in  x ••  27.86  (in  x)2 

160 

y = 280.51  - 141.90  in  x + 23.96  (in  x) 2 

-65 

y = 544.94  - 210.69  in  x+  21.69  (in  x) 2 

2 

70 

y = 333.50  - 150.03  in  x + 18.71  (in  x) 2 

160 

y = 202.97  - 93.58  in  x+  13.17  (in  x) 2 

-65 

y = 517.16  - 194.42  in  x + 18.94  (in  x)“ 

3 

70 

y = 289.25  - 123.01  in  x + 13.98  (in  x) 2 

160 

y = 170.46  - 73.58  in  x + 9.21  (in  x)2 

Appendix  B 

STEPWISE  REGRESSION  PROGRAM  LISTING 
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Figure  C-l.  -65°F,  12-inch  drop  height 


Figure  02.  -65°F,  24-inch  drop  height 
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Figure  C-5.  160°F,  12-inch  drop  height 
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Figure  D-l.  -65°F,  18-inch  drop  height 
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Figure  D-3.  70°F,  18-inch  drop  height 


Figure  D-4.  70°F,  30-inch  drop  height 
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Figure  D-5.  160°F,  18-inch  drop  height 
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Figure  D-6.  160° F,  30-inch  drop  height 
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Figure  E-l.  -65°F,  12-inch  drop  height 


Figure  E-2.  -65°F,  24-inch  drop  height 
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Figure  E-4.  70°F,  24-inch  drop  height 


Figure  E-5.  160°F,  12-inch  drop  height 


Figure  E-6.  160°F,  24-inch  drop  height 
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TWENTY- SEVEN-INCH  DROP  HEIGHT  MINICEL  DATA, 
1,  2,  AND  3 INCH  THICK 


132 


m 


0.20  0.80 


Appendix  G 

POUR-  AND  FIVE -INCH  THICK  MINI CEL  DATA 
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MINICEL  - 18  in.  Drop  Height  (Continued) 
STRESS  LEVELS  (PSI) 
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STRESS  LEVELS  (PSI) 
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MINICEL  - 24  in.  Drop  Height  (Continue 
STRESS  LEVELS  (PS1) 


o 

C*) 

r-H 

CNI 

CNI 

CNI 

in 

m 

CNI 

CNI 

CM 

cn 

O 

o> 

MINICEL  - 30  in.  Drop  Height  (Continued) 
STRESS  LEVELS  (PSI) 


h r j n h cni  n 

U0T^H0T-[d3^ 


Appendix  H 

PROGRAM  LISTING  OF  CUSHION  OPT 
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HAS  THIRO  TEMPERATURE  BEEN  REACHEO 
N .GE.  3 > GO  TO  500 
INCREMENT  TEMPERATURE 
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